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External driving is emerging as a promising tool for exploring new phases in quantum systems. 

The intrinsically non-equilibrium states that result, however, are challenging to describe and control. 

We study the steady states of a periodically driven one-dimensional electronic system, including the 
effects of radiative recombination, electron-phonon interactions, and the coupling to an external 
fermionic reservoir. Using a kinetic equation for the populations of the Floquet eigenstates, we 
show that the steady state distribution can be controlled using the momentum and energy relaxation 
pathways provided by the coupling to phonon and Fermi reservoirs. In order to utilize the latter, 
we propose to couple the system and reservoir via an energy filter which suppresses photon-assisted 
tunneling. Importantly, coupling to these reservoirs yields a steady state resembling a band insulator 
in the Floquet basis. The system exhibits incompressible behavior, while hosting a small density of 
excitations. We discuss transport signatures, and describe the regimes where insulating behavior is 
obtained. Our results give promise for realizing Floquet topological insulators. 


The availability of coherent driving fields such as lasers 
opens many exciting possibilities for controlling quantum 
systems. In particular, the recent realization that the 
topological characteristics of Bloch bands can be modi¬ 
fied through periodic driving® sparked a wave of pro¬ 
posals^ and experiments^^ to realize various types 
of “Floquet topological insulators” in solid state, atomic, 
and photonic systems. Here topology emerges in the 
basis of Floquet states, time-periodic eigenstates of the 
driven system’s single-particle evolution operatoi^HUl 

Floquet states provide a convenient basis for describing 
the evolution of driven systems, in many ways analogous 
to the Hamiltonian eigenstates of non-driven systems. 
However, the powerful thermodynamic rules that govern 
the level occupations of static systems in thermal equi¬ 
librium in general cannot be directly translated into the 
inherently non-equilibrium context where Floquet states 
are defined®®^. Photon-assisted scattering processes in 
which energy is exchanged with the driving field produce 
heating and violate the reversibility conditions that give 
rise to the Boltzmann distribution in equilibrium 27 . This 
crucial difference brings up many intriguing and impor¬ 
tant questions about the steady-state physical properties 
of open Floquet-Bloch systems. In particular, in order to 
realize the promise of non-equilibrium topological phe¬ 
nomena, one of the major outstanding problems is to 
identify which types of systems, baths, and system-bath 
couplings can lead to non-equilibrium steady states en¬ 
abling Floquet topological insulators to exhibit b ehavio rs 
similar to those of their equilibrium counterpart^ - 

Floquet-Bloch steady-state distributions are currently 
known for certain special cases. A closed driven system 
inevitably heats up t o a m aximal entropy (infinite tem¬ 
perature) steady stat d®^ , in the absence of many-body 
localizatioi l — I—I In contrast, an open driven system con¬ 
nected to a thermal bath need not reach such an end. In¬ 
deed, when the system and system-bath coupling Hamil¬ 
tonians can be made time-independent through a simple 
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FIG. 1. Carrier kinetics in a Floquet-Bloch system coupled to 
Bose and Fermi reservoirs, a) One dimensional semiconductor 
wire coupled to an energy-filtered fermionic reservoir. Energy 
filtering is achieved by coupling the system and reservoir via 
a deep impurity band in a large bandgap semiconductor, b) 
Band structure of the non-driven system. The driving field 
photon energy Ml exceeds the bandgap -E gap , causing reso¬ 
nant coupling at crystal momentum values ±fci?. c) Floquet 
band structure, indicating the character of the Floquet band 
in terms of the original conduction (blue) and valence (red) 
bands. Coupling to acoustic phonons mediates electronic mo¬ 
mentum and energy relaxation (orange arrows), while radia¬ 
tive recombination scatters electrons vertically between con¬ 
duction and valence band like states (purple arrow). At half 
filling, the steady state resembles that of an insulator with a 
small density of excited electrons and holes. 

rotating frame transformation, Gibbs-type steady states 
are expectecP®^. More generically, however, even spon¬ 
taneous emission into a zero temperature bath may cause 
heating due to the possibility of absorbing energy from 
the driving field. How to control the steady states of 
driven systems, and, in particular, under what conditions 
they may be used to explore novel topological phenom- 
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ena, remain a challenging problem. 

Our aim in this work is to uncover new means to con¬ 
trol the steady state occupations of Floquet-Bloch states 
in driven systems. Here we focus on the dissipative 
open-system dynamics governed by the system’s coupling 
to external baths; the complicated problem of electron- 
electron interactions will be addressed in future work. 
In particular we target the case of a half-filled fermionic 
system, where we seek to obtain an insulator-like steady 
state in which the lower Floquet band is filled and the 
upper Floquet band is empty. We refer to this state as a 
Floquet insulator. We investigate how this state can be 
approached through the relaxation of momentum and en¬ 
ergy, enabled by connection to low temperature bosonic 
and fermionic baths, see Fig. [T] For a semiconductor- 
based realization, these baths naturally correspond to 
phonons and the electromagnetic environment (bosonic 
baths), and to a lead connected to the system (fermionic 
bath); analogo us co uplings can be arranged, e.g., in cold 

atomic system^S122] 


Several dynamical processes and their interplay gov¬ 
ern the density of particle-hole “excitations” above the 
ideal Floquet insulator stat^I. Radiative recombination 
constitutes an important mechanism for generating ex¬ 
citations. In a non-driven system, recombination allows 
electrons in the conduction band to annihilate with holes 
in the valence band via the spontaneous emission of a 
photon. For resonant driving, as illustrated for the case 
of a one-dimensional system in Fig. Eh the Floquet bands 
feature a band inversion: states with crystal momenta be¬ 
tween the two resonance values ztkp in the lower Floquet 
band are primarily formed from conduction band states 
of the non-driven system, while in the same interval the 
upper Floquet band is comprised of valence band states. 
Therefore, radiative recombination results in transitions 
from the lower to the upper Floquet band, thus increas¬ 
ing the density of excitations, see Fig. [jj:. Phonon scat¬ 
tering, on the other hand, enables relaxation of momen¬ 
tum and quasi-energy within and between Floquet bands, 
and may balance the recombination-induced heating. A 
fermionic reservoir provides additional channels for re¬ 
moving excitations from the system, and also gives means 
to tune its total carrier density. Importantly, photon- 
assisted electron-phonon scattering, as well as photon- 
assisted tunneling to/fro m the Fermi reservoir, generally 
also contribute to heatin^fiMU, see Fig. [2j 


Our main message is that the driven electronic system 
can approach the Floquet insulator steady state when 
appropriately coupled to phonon and Fermi reservoirs. 
In order for this to work, the coupling to the fermionic 
reservoir must be “engineered” to avoid the deleterious 
effects of photon-assisted tunneling. This can be accom¬ 
plished by connecting the system to the reservoir via a 
narrow-band energy filter (see Fig. [2] and Sec. IIIC). We 
also discuss regimes in which photon-assisted electron- 
phonon scattering can be suppressed. Most remarkably, 
at low temperatures and with energy-filtered coupling 
to a fermionic reservoir, we find that the driven system 


exhibits incompressible and insulating behavior. This 
implies that a steady-state Floquet topological insulator 
phase may be within reach. 


A. Structure of the paper and main results 

Before beginning the analysis, we briefly summarize 
the structure of the text to follow. Keeping in mind 
our motivation of realizing Floquet topological insula¬ 
tors, our main focus in this work is on achieving Floquet 
insulator steady states. 

First, in Sec. [I] we introduce the Floquet states of the 
periodically-driven lattice system, with band structure 
depicted schematically in Fig. |T|>. After defining the 
Floquet states, we introduce the Floquet kinetic equa¬ 
tion, Eq. 0, which forms the basis for the description of 
many-body population dynamics throughout this work. 
The kinetic equation can be obtained systematically from 
the exact (infinite) hierarchy of equations of motion for 
multi-particle correlators (see Appendix [A|, and at our 
level of approximation takes on a simple intuitive form 
in terms of incoming and outgoing rates for each state. 

Next, in Sec. |TT] we study the steady states when the 
system is coupled only to the bosonic baths. Here our 
aim is to elucidate the competition between heating due 
to radiative recombination and momentum/energy relax¬ 
ation by phonons in a particle number conserving sys¬ 
tem. When the electron-phonon scattering rates (ignor¬ 
ing Pauli blocking) are large compared to the recombi¬ 
nation rate, we find that the driven system approaches a 
Floquet insulator state, with separate particle and hole 
densities in the upper and lower Floquet bands, respec¬ 
tively, see Fig. [3] below. The steady state excitation 
density depends on the ratio of phonon-assisted inter- 
Floquet-band relaxation and recombination rates, be¬ 
coming small for fast interband relaxation. Using rather 
general arguments, we show that the steady state exci¬ 
tation density scales with the square root of the recom¬ 
bination rate in the limit of fast interband relaxation. 
As a result, even strong electron-phonon coupling may 
be insufficient to fully deplete excited carriers from the 
system. 

In Sec. |III| we introduce coupling to a fermionic reser¬ 
voir. In Fig. [4] we display the steady states for both 
wide-band and energy-filtered reservoirs. We show that 
coupling to a wide-band reservoir increases the density of 
excitations, due photon-assisted tunneling. The energy- 
filtered system-reservoir coupling suppresses all photon- 
assisted tunneling, and our results demonstrate that it 
can further reduce the density of excitations. We discuss 
two coupling geometries, where the Fermi reservoir is ei¬ 
ther coupled to the system at a single point (as a lead), or 
uniformly throughout the system. For homogeneous cou¬ 
pling, when the chemical potential of the filtered reser¬ 
voir is set inside the Floquet gap, the excitation density 
may be highly suppressed, thus bringing the system close 
to the ideal Floquet insulator state. Interestingly, even 
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when the steady state hosts a finite density of excitations, 
the system is incompressible in the sense that the excita¬ 
tion density is unaffected by small shifts of the chemical 
potential of the reservoir, see Fig. [5] For a lead coupled 
at a single point, the steady state distribution is neces¬ 
sarily inhomogeneous. We provide an estimate for the 
“healing length” over which the distribution can be af¬ 
fected by such coupling. Beyond this length, the steady 
state is set by the competition between recombination 
and electron-phonon coupling, as described in Sec. [II] 
Finally, in Sec. IV we summarize the main results and 
discuss implications for transport experiments. We dis¬ 
cuss the corresponding observables and the conditions 
under which insulating behavior could be observed. 


I. FLOQUET-BLOCH KINETIC EQUATION 
FOR THE DRIVEN TWO-BAND SYSTEM 

In this section we describe the single-particle proper¬ 
ties of an isolated periodically-driven system. We first 
give the Hamiltonian of the system without driving, and 
then discuss the form of driving and the resulting Flo- 
quet states. We then introduce the kinetic equation for 
Floquet-state occupation numbers, which is the foun¬ 
dation for the description of nrany-body dynamics used 
throughout this work. The section concludes with a brief 
overview of the dynamical processes described by the ki¬ 
netic equation. 


A. System Hamiltonian and Floquet-Bloch states 

We now introduce the single-particle Hamiltonian and 
Floquet-Bloch states for the periodically-driven system 
that we consider. Many of the features that we describe, 
including the form of the kinetic equation, hold quite gen¬ 
erally, independent of dimensionality. For concreteness, 
and to allow comparison with detailed numerical simula¬ 
tions, we focus on the case of a one-dimensional system 
with two bands. 

The single particle Hamiltonian of the driven system 
is defined as follows. We assume that the driving field 
is spatially uniform, thus maintaining the translational 
symmetry of the lattice. In this case the crystal momen¬ 
tum k is conserved. For each k the evolution within the 
corresponding 2x2 Bloch space is given by the Hamil¬ 
tonian H(k) = H 0 (k ) + V(t), with 

H 0 (k) = \E k {d k ■ a), V(t) = \Vq {g ■ <t) cosfli, (1) 

where i^E k are the energies of the conduction and va¬ 
lence bands, d k and g are unit vectors, Vo and Q are the 
amplitude and angular frequency of the drive, and <x is 
a vector of Pauli matrices acting in the two-component 
orbital space (in this work we ignore spin). For now we 
leave the values of d k and g unspecified, giving explicit 
forms when discussing numerical results below. 


To understand the nature of the coupling induced by 
driving, we rotate to the basis of conduction and va¬ 
lence band states, i.e., to the basis which diagonalizes 
H 0 (k). Specifically, the Bloch eigenstates in the conduc¬ 
tion and valence bands of the non-driven system satisfy 
dlo(^)l^fcc) — Ho{E)\v, kv ) — 2 | ^kv ) • 

The driving term V (t) in Eq. 0 is expressed in the ba¬ 
sis of lattice orbitals, and naturally does not depend on 
k. However, after rotating to the basis of conduction and 
valence band states for each k, the driving picks up a 
non-trivial ^-dependent matrix structure 

H 0 (k) = \E k a z , V(k,t) = \V 0 (g k • er) cos fit, (2) 

where tildes indicate operators in the basis of conduction 
and valence bands, and g k = g k .\\z + 9k,± is a unit vector 
determined by the relative orientations of d k and g in 
Eq. 0 . broken down to z and x — y components. 

When the system is isolated, the Floquet-Bloch states 
{\ipk±(t))} provide a convenient basis for describing its 
evolution. Each state \ipk±(t)) can be expressed as a sum 
over harmonics: 

OO 

|^±(*))= E (3) 

n =—oo 

where £ k ± is the quasi-energy of | i/) k ±(t)) and {\(j) k± )} 
is a non-normalized (and over-complete) set of states 
found by Fourier transforming the time-dependent 2x2 
Schrodinger equatiorPfil- 1 in the Bloch space for crystal 
momentum k. The quasi-energies {£ k ±} and harmon¬ 
ics {|in Eq. (|]) are only uniquely defined up to 
the gauge freedom £' k ± = £ k ± + mMl, \<j/£±) = |< i>^± n )- 
Here we fix the gauge by choosing £ k ± within a single 
Floquet-Brillouin zone centered around a specific energy 
Eq, Eq — Ml/2 < £ k ± < Eq T Ml/2. 

Before discussing many-body dynamics, a few com¬ 
ments on the nature of the Floquet bands are in order. 
We are interested in the case where the driving field pho¬ 
ton energy Ml exceeds the bandgap -E gap of the non- 
driven system, see Fig. [I]>. In the Floquet picture, the 
leading-order influence of the driving can be understood 
by first shifting the valence band up by the photon en¬ 
ergy Ml. After shifting, the bands become degenerate at 
the resonance points 42 ±kf j in the Brillouin zone where 
E kR = Ml. Here the driving opens avoided crossings, 
resulting in a gap A kR ss Vo|9fc H ,_L| between the two Flo¬ 
quet bands. The resulting band structure is depicted in 
Fig- [TJ • We center the Floquet zone on these resonances 
in the conduction band, setting E 0 = \Ml. Throughout 
we assume that the bandwidth is narrow enough such 
that the two-photon resonance condition is never satis- 
fiecP®, i.e., E k < 2hfl for all k. 

As discussed in the introduction, the resonant driving 
introduces a band inversion in the Floquet bands. Fur¬ 
thermore, near the resonant momenta kn the Floquet 
bands are strongly hybridized superpositions of conduc¬ 
tion and valence band states. These features of the Flo¬ 
quet bands have important consequences both for con- 
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trolling band topolog}®^ and for the nature of many- 
body dynamics in the system, as we describe below. 


B. The Floquet kinetic equation 

Below we use the Floquet basis of single-particle states 
to describe the many-body dynamics of the driven system 
when it is coupled to bosonic and fermionic baths. In 
particular, we aim to characterize the steady states of the 
system in terms of the Floquet state occupation numbers 
F ka = (ft a (t)fka{t)) 7 where f ka {t) creates an electron in 
the state | if>ka(t)) at time t, with a = ±. Focusing on 
the dynamics for time scales much longer than the driving 
period, we develop a kinetic equation in the Floquet basis 
to describe the net rate of change of the population in the 
Floquet state | ipka) due to electron-phonon scattering, 
radiative recombination, and tunneling to and from the 
fermionic reservoir: 

Fka = i£({F}) + r£({F}) + It™ (F ka ). (4) 

Here the “collision integrals” / ph , 7 rec , and J tun describe 
electron-phonon scattering, recombination, and tunnel 
coupling to the reservoir, respectively, and {F} indicates 
the set of occupation factors for all momentum and band 
index values. The key processes associated with each of 
these terms are represented schematically in Fig. [l}\ 

The derivation of Eq. ([4]) is rather technical, so here 
we briefly summarize the approach (for details, see Ap¬ 
pendix [A] and, e.g., Ref. [46j) . We begin by consider¬ 
ing the equations of motion for the single-particle cor¬ 
relators ■ Coupling to the bath degrees of 

freedom generates an infinite hierarchy of equations of 
motion involving correlators of higher and higher order. 
We focus on a homogeneous system, in the regime where 
coherences between different Floquet states can be ne¬ 
glected (see below). Using a standard cluster-expansion 
approach, we systematically truncate the equation of mo¬ 
tion hierarchy and obtain transition rates which coincide 
with those given by the “Floquet Fermi’s golden rule.” 
Below we frame the discussion in terms of these golden- 
rule transition rates, which provide a clear intuitive pic¬ 
ture for all terms contributing to Eq. Q. We will use 
the rates to build up the specific forms of the collision 
integrals / ph , J rec , and J tun . 

In describing the dynamics of the system, it is impor¬ 
tant to note that the occupation factors F^a do not gener¬ 
ally give a complete description of the steady state. How¬ 
ever, when transition rates associated with the system- 
bath interaction are small in comparison with the Floquet 
gap off-diagonal correlations such as (ft„(t)fk a '(t)) 
are suppressed in the steady state (see Appendix[A I. Cru¬ 
cially, even if the scattering rates are large when the sys¬ 
tem is far from the steady state, Pauli blocking in the 
steady state can strongly suppress the phase space for 
scattering. Therefore Eq. © can provide a good descrip¬ 
tion of the steady state, even in parameter regimes where 


it does not give a faithful description of the transient dy¬ 
namics. The requirement that the steady state scattering 
rates remain small compared with the Floquet gap 
provides an important consistency check, which we ap¬ 
ply to all numerical simulations discussed below. Note 
that the cluster expansion approach provides a powerful 
framework that can be used to incorporate the roles of 
coherences and non-Markovian dynamics, going beyond 
the regime studied here 47 . 


II. ELECTRON-PHONON COUPLING AND 
RECOMBINATION 

In this section we discuss the steady states of the elec¬ 
tronic system which result from the competition between 
radiative recombination and coupling to the phonon 
bath. Both processes arise from the coupling of elec¬ 
trons to a bosonic bath, comprised of photons in the 
former case and phonons in the latter. Formally, the 
collision integrals J rec and I ph describing these processes 
are very similar. However, it is important to understand 
that they act in competition. During recombination, an 
electron transitions from the non-driven conduction band 
to the valence band, while emitting a photon. In terms 
of the Floquet bands, this process promotes an electron 
from the lower to the upper Floquet band (see Fig. [l]), 
thereby heating the electronic system and increasing the 
total number of excitations. On the other hand, electron- 
phonon scattering primarily relaxes excited electrons to 
the bottom of the upper Floquet band (and similarly re¬ 
laxes holes to the top of the lower Floquet band), and 
allows excited electrons to relax back to the lower Flo¬ 
quet band, thereby reducing the number of excitations. 

Note that the electron-phonon interaction may also 
play an adverse role in the system: photon-assisted scat¬ 
tering processes may increase the number of excitations. 
We show that these processes can be effectively elimi¬ 
nated under suitable conditions on the phononic spec¬ 
trum and the form of the drive. Even when these pro¬ 
cesses are eliminated, radiative recombination remains as 
a source of heating in our model. 

The competition between electron-phonon scattering 
and recombination determines the steady state of the 
system. We show that these steady states feature Fermi 
seas of excited electrons and holes, with separate chem¬ 
ical potentials and a temperature equal to that of the 
phonon bath. Below we first analyze the kinetic equa¬ 
tion in the presence of a generic bosonic bath. We then 
input the specific details needed to describe recombina¬ 
tion and scattering by acoustic phonons, and analyze the 
resulting steady states. Finally, through an analytical es¬ 
timate we show that for fixed electron-phonon coupling 
the steady state excitation density grows with the square 
root of the radiative recombination rate. 





5 


A. Collision integral for a generic bosonic bath 


The bosonic bath is described by the Hamiltonian 
Hb = Y^qhuq^qbq, where b' q and b q are the creation 
and annihilation operators for a bosonic excitation car¬ 
rying (crystal) momentum q , and ui q is the corresponding 
frequency. Using the creation and annihilation operators 
{c(. v , c ku} for electrons in the bands of the non-driven 
system, defined below Eq. |l]) , we describe the “electron- 
boson” interaction via H\ ni = Y q H- ln t(q), with 

tfint(g) = E E G t"\<l)ciyCUb\ + b_q). (5) 

kk' ip,iy' 

Here G\” (q) is the matrix element for scattering an elec¬ 
tron with crystal momentum k in band v to crystal mo¬ 
mentum k' in band v', with the emission (absorption) 
of a boson of momentum q (—q). We take the bath to 
be three dimensional, and the electronic system to lie 
along the x-axis. Note that in Eq. (|5j) we did not im¬ 
pose k' = k — q x to allow the possibility of describing 
a finite system coupled to a bath of larger size (such as 
the electromagnetic environment). If the lengths of the 
bath and system are the same, we can impose conserva¬ 
tion of the corresponding crystal momentum component, 
whereby G\ff (q) is nonzero only if k! = k — q x . 

As a fundamental building block for constructing the 
many-body collision rates, we calculate the rate for 

a single electron in an otherwise empty system to scatter 
from crystal momentum k in Floquet band a to crystal 
momentum k' in Floquet band a'. For transparency, we 
focus on zero temperature; the analogous expressions at 
finite temperature are given in Appendix [A] For simplic¬ 
ity we take the matrix elements in Eq. |5]) to depend only 
on q x , i.e., G\ff (q) = G\” (q x ); the discussion that fol¬ 
lows can be easily generalized beyond this assumption, 
but the qualitative results will not be affected. 

Due to the harmonic structure of the time-dependent 
Floquet state wave functions, Eq. ([3]), the transition rate 
is given by a sum over many contributions, W= 
(n). In terms of the electronic operator G(q x ) = 

Sfc,fc' Yj V y G ku '(. ( lx) c k', I s c kv'i these contributions are 
given by 


Wf{n) 


?7T I 

Y E\E^\G(q x MZ) 

q x rn 


PqA~^ £ n), 


( 6 ) 

where A£ n = £k' a '—£ka + nhCl. The quasi-energy 
difference between final and initial electronic states is 
£k'a'—£kon and PqS^ 1 ) is the boson density of states at 
frequency u> for a fixed value of the boson’s longitudinal 
momentum component q x . Note that for a monotonic bo¬ 
son dispersion, p 9x (w) is only nonzero if u) > w qo , where 
Qo = (qx, 0,0). The scaling of the individual rates Wjfff 
with system size is discussed in Appendix |B| 

The structure of the transition rates in Eq. © can 
be understood heuristically as follows. Superficially, 
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FIG. 2. Harmonic structure of Floquet states and energy- 
filtered reservoir coupling, a) Floquet harmonics of a two-level 
system with states |v) and |c) coupled by an on-resonance 
driving field V(t). The Floquet zone (shaded) is centered at 
the energy Eq, set equal to the energy of the resonant state |c). 
In the special case of a rotating-field, V ( t) = | VQe~' nt \c){v\ + 
h.c., we have |0±) = |c), I^J 1 ) = ±|u) and £± = Eq±^Vo, see 
Eq. ©• Away from resonance, the relative normalizations of 
\<f>+) and \cf)-) will change. For a more general form of weak 
driving, the dominant harmonics are shown in bold, b) The 
Floquet states |^>±(£)) are both coupled to filled and empty 
state s of a wide-band reservoir via the harmonics {| <^>±)}, see 
Eq. (101. Here the reservoir chemical potential is set in the 


gap of the non-driven system, c) When coupling is mediated 
by a narrow-band energy filter, the tunneling density of states 
(TDOS) and photon-assisted tunneling are suppressed outside 
the filter window. By setting the reservoir chemical potential 
inside the Floquet gap, centered around the energy Eq in the 
original conduction band (see Fig. the lower and upper 
Floquet bands are selectively filled and emptied, respectively. 


Hk±(t)) in Eq. © takes the form of a superposition over 
a ladder of states \(K±) with “energies” £k± + nKl, see 
illustration in Fig. [2a. Viewing these harmonics \4>k±) 
as independent states, the net transition rate = 

Y^ n ^ka C n ) f° unc l by summing the contributions from 
all pairs of initial and final states, while taking into ac¬ 
count “energy” conservation. The appearance of nhfl 
inside the density of states in Eq. © expresses the fact 
that quasi-energy is a periodic variable, and therefore 
Floquet scattering processes need only conserve quasi¬ 
energy up to multiples of the driving field photon energy 
hQ. If a boson is emitted and an electron decreases its 
quasi-energy, A£ < 0, then the scattering rate can be 
non-zero for n = 0. Interestingly, the scattering rate can 
also be non-zero if a boson is emitted and an electron 
increases its quasi-energy, AS > 0, if n < 0. 

The collision integrals in Eq. Q are given by the differ¬ 
ences between the total rates for scattering into and out 
of the state \ifka(t)), due to recombination or coupling 
to acoustic phonons. In turn, these rates are obtained by 
multiplying the bare rates in Eq. ([6]) by products of oc¬ 
cupation factors Fka, Fk'a’ = (1 — Fk> a '), etc., to account 
for the filling of the initial and final states: 


h a = E {w^FkaFk'a' - Wka Fk'a’F k<x 

k' ol' 


( 7 ) 


The corresponding expressions for nonzero bath temper- 
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ature are shown in the Appendix [A] 

Equations Q and ([7|) support what we refer to as 
“Floquet-Umklapp” processes, in which quasi-energy 
conservation is satisfied with n ^ 0. Such processes 
generically heat the system when they are allowed within 
the kinematic constraints imposed by the bath and 
Floquet-system spectra (i.e. by energy and momentum 
conservation). For example, even at zero bath tempera¬ 
ture, an electron may spontaneously jump from the lower 
Floquet band to the upper one while emitting a bosonic 
excitation (see Fig. So- As we show below, such pro¬ 
cesses cause deviations from the ideal Floquet insulator 
state. 

Fortunately, Floquet-Umklapp processes are sup¬ 
pressed under appropriate conditions on the dispersion 
of the bath bosons. In fact, Floquet-Umklapp processes 
are completely suppressed if the bath bandwidth is lim¬ 
ited such that Wjfff ( n ) strictly vanishes for all n / 0. 
Practically speaking, this means that the maximal bo¬ 
son energy must be smaller than the quasi-energy 

gap at the Floquet zone edge (i.e., the gap around Ml/2 
in Fig. such that the energy conservation condition 
—AS — hu .;™ ax — nMl = 0 cannot be satisfied with « / 0. 
Below we will show how Floquet-Umklapp processes are 
manifested in radiative recombination and phonon scat¬ 
tering processes, and discuss methods to suppress them. 


B. Radiative recombination 

Having established the general framework for coupling 
the driven system to a bosonic bath, we now use it to 
study specific dissipation mechanisms which are relevant 
for driven semiconductor systems. We start by consider¬ 
ing radiative recombination. 

In non-driven systems, radiative recombination occurs 
when an excited particle in the conduction band relaxes 
to fill a hole in the valence band. This results from the 
interaction of electrons with the electromagnetic envi¬ 
ronment, which is represented by a bosonic bath in our 
model. In typical semiconductors, the electromagnetic 
interaction leads to transitions between states of differ¬ 
ent bands. This restriction on the transitions arises due 
to two important factors: 1) the large speed of light im¬ 
plies that energy and momentum conserving transitions 
are practically “vertical” (i.e. the electronic momentum 
is conserved), and 2) the electromagnetic dipole matrix 
elements couple states from different atomic orbitals. To 
impose this restriction in our model we describe the inter¬ 
action with the electromagnetic environment using ma¬ 
trix elements of the form G rec oc (1 — £„„>). For simplicity, 
in the simulations below we model vertical recombination 
vigP0G rec = ff rec (l — S vv i)6 qmfi 6 k>k >, and take a constant 
density of states po for photons with energies Hco > E gap . 

We now describe the processes resulting from the cou¬ 
pling to the electromagnetic environment in the driven 
system that we consider. Following from the situation in 
the non-driven case, relaxation via emission of a photon 


to the environment is possible from a Floquet state of 
predominantly conduction band character to one of pre¬ 
dominantly valence-band character. Due to the band in¬ 
version described in detail in Sec. |TJ the — Floquet band 
has conduction band character for momenta fc < k r. 
Therefore, spontaneous transitions from the — to the + 
Floquet band are possible for states in this momentum 
range. Note that these Floquet-Umklapp processes in¬ 
crease the total electronic quasi-energy, and play an im¬ 
portant role in determining the density of excitations in 
the steady state of the system (see Sec. IID). The rates 
of these processes may be controlled to some extent by 
placing the system in a cavity or photonic crystal, which 
modifies the photon density of states. In addition, spon¬ 
taneous transitions from the + to the — Floquet band are 
allowed in the momentum region \k\ > Ur, where the + 
Floquet band has predominantly conduction band char¬ 
acter. These processes help to reduce the total electronic 
quasi-energy, but will play a minor role near the steady 
state where the + Floquet band is mostly empty. 

The processes we have considered so far follow directly 
from those that are active in a non-driven system. How¬ 
ever, in a driven system an electron may also transition 
from a state of valence band character to one of conduc¬ 
tion band character, by emitting a photon to the envi¬ 
ronment while absorbing energy from the drive. Such 
processes are possible for initial states in the — Floquet 
band with |/c| > kf j, and for initial states in the + Floquet 
band with fc < /cr. Referring to Eq. Q, the matrix ele¬ 
ments for these processed are suppressed by [Vo /{hVL)} 2 
for weak driving, and hence their rates are suppressed as 

[v 0 /(Mi))\ 


C. Scattering due to acoustic phonons 

The interaction between the electronic system and a 
bath of acoustic phonons plays a key role in setting the 
steady state of the driven system. Phonon-mediated scat¬ 
tering quickly relaxes excited electrons (holes) to the bot¬ 
tom (top) of the respective Floquet band. In addition, 
phonon-mediated scattering allows these excitations to 
relax across the Floquet gap. The competition between 
the latter interband scattering processes and radiative re¬ 
combination sets the steady state density of excitations, 
as we discuss below. 

In our model we assume that the electron-phonon 
coupling conserves the band index v of the non-driven 
system, G ph oc 5 VI J >, as is typical for wide gap 
semiconductor jll ?l 50 l 51 l The coherent drive hybridizes the 
bands near the resonances ±/cr, thus enabling both in¬ 
traband and interband scattering in the Floquet bands 
(see Fig. Eo- Note that the scattering crucially involves 
the exchange of both crystal momentum and quasi¬ 
energy between the phonons and the electrons, thus al¬ 
lowing relaxation of these quantities. We take the ma¬ 
trix elements to conserve lattice momentum, G\f ( q x ) = 
g(ijx)d vv 'S, h:k -k>■ In principle, the q x dependence of 
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g{q x ) depends on the specific type of electron-phonon 
coupling. For simplicity, we take the matrix elements 
to be independent of q x , but have numerically verified 
that other choices do not change the qualitative results. 

Besides helping to relax excitations, photon-assisted 
electron-phonon scattering can increase the excitation 
density. Such Floquet-Umklapp scattering transfers elec¬ 
trons from the lower to the upper Floquet band, and can 
occur even for a zero temperature phonon bath. 

Phonon-related Floquet-Umklapp processes can be 
suppressed in several ways. First, as discussed in 
Sec. m limiting the bandwidth for the phonons to 
be smaller than the quasi-energy gap at the Floquet 
zone edge, A ec j ge , efficiently suppresses photon-assisted 
scattering. Note, however, that the phonon bandwidth 
should remain bigger than the Floquet gap emerging 
at the resonance momenta, Afc fl , as otherwise phonons 
would be unable to facilitate relaxation between the up¬ 
per and lower Floquet bands. An optimal phonon band¬ 
width u>d would therefore satisfy Ak R < Hlud < A 0 d ge - 
The bandwidth for the phonon bath depends on material 
parameters, however, and may not be easily tunable. 

Interestingly, additional routes are available for sup¬ 
pressing Floquet-Umklapp processes involving phonons. 
If the boson bandwidth allows the energy conservation 
condition — AS — Sw“ ax — nhfl = 0 to be satisfied for 

\n\ < 1 (but not for \n\ > 1), the rates W^(n) with 
n = ±1 can be controlled by the choice of driving. In 
particular, for harmonic driving they vanish as gy 0 
(see Appendix[A]). For many experimentally-relevant ma¬ 
terials driven by optical fields, g\\ is indeed small for mo¬ 
menta near k = 0. Additionally, even when none of the 
conditions above are met, the amplitudes of the Floquet 
harmonics {| </>£„)} (and hence the rates (n)) are 
generically suppressed for large n. Hence, although heat¬ 
ing inevitably accompanies coupling to a bosonic bath, 
there are many ways to control or limit the correspond¬ 
ing effects on the steady state distribution (see below and 
also Refs. I53H57|) . 


D. Steady state 


The steady state of the driven model described above 
results from the competition between the two main dissi¬ 
pation mechanisms: radiative recombination and acous¬ 
tic phonon scattering. To gain a more quantitative pic¬ 
ture of the behavior, we numerically solve for the steady 
states of the kinetic equation @ in the model out¬ 
lined above, with the parameter values given in Table [Tj 
We take acoustic phonons to have a linear dispersion 
in three dimensions, CtT h = c s |q|, up to a “Debye fre¬ 
quency” cutoff W£>. We focus on the situation where 
Afc H < Hlod < A ec jge, which allows inter-Floquet-band 
scattering, but forbids Floquet-Umklapp phonon scat¬ 
tering processes. The rates {W^,} are calculated us¬ 
ing the form for the matrix elements described in Sec¬ 
tions IIB and IIC Our results are summarized in Fig. [3] 


A 

E gap 

g 

d k 

Vo 

o.25 hn 

o.8 nn 

(1,0,0) 

(0,0,1) 

0.1 Ml 

Cs 

huJD 

2tt (G p 0 h ) 2 p ph 

kgT 

7T \/3 

0.15/iD 

(2 x lO" 2 ) hQ 

0.1Afc H 


TABLE I. Parameters fixed in all simulations. Top row: 
parameters of the electronic Hamiltonian, Eq. (JT]) , with E = 
2A[1 — cos(fca)] + E gap , where a is the lattice constant. The 
drive is spatially uniform, V(t) = dUo(g ■ cr) cos Of. Bottom 
row: parameters of the three dimensional acoustic phonon 
bath, where c s is the phonon velocity, and wn is the Debye 
frequency. In all simulations, the overall scale of the phonon 
matrix elements is set by fixing the ratio 2n(GQ h ) 2 p ph /(fifl), 
where p ph is the phonon density of states at zero momentum 
and energy hc 3 {iY/a). For convergence, in the simulations we 
keep the phonon bath at a small temperature, fesT « 10 _2 ftfh 

The main result of the numerical simulation is the 
steady state distribution of excited carriers in the 
upper Floquet band, which is shown in Fig. [3]a,, for a to¬ 
tal density of particles corresponding to half-filling. Due 
to particle-hole symmetry, the distribution of holes in 
the lower Floquet band, Ffc_, is identical to that above. 
We examine the behavior of the steady state distribu¬ 
tion while tuning the ratio of the phonon scattering and 
radiative recombination rates. Specifically, we fix the 
parameters for electron-phonon coupling, and vary the 
overall scale of the matrix elements for recombination. 

As seen in Fig. [3}i. in all cases the upper Floquet band 
hosts a finite density of excited Floquet carriers, local¬ 
ized around the two band minima. For relatively weak 
electron-phonon coupling the excitation density is large, 
but is limited by saturation of the recombination rates 
due to Pauli exclusion above a given density. Notably, 
when electron-phonon coupling is relatively strong, the 
excitation density is significantly suppressed. Moreover, 
in this regime the distribution of excited carriers is well 
described by a Floquet-Fermi-Dirac distribution with an 
effective chemical potential p, e (a fit parameter) and tem¬ 
perature corresponding to that of the phonon bath (solid 
lines). By “Floquet-Fermi-Dirac distribution” we refer to 
a distribution of particles in Floquet states, which is de¬ 
scribed by a Fermi-Dirac distribution taken as a function 
of quasi-energy. The distribution of holes in the lower 
Floquet band (not shown) takes an identical form with 
an effective chemical potential /.% = p e due to particle- 
hole symmetry of the model. To check the consistency of 
our approach, we verify that the scattering rates in the 
steady state are significantly smaller than the Floquet 
gap Afc B . This condition is satisfied in particular for mo¬ 
menta around kn where the electron and hole excitation 
densities are localized. A more detailed discussion of the 
scattering rates is provided in Appendix |D| 

The above form for Fk + can be understood by consid¬ 
ering the dynamics of the electrons coupled to the pho¬ 
ton and phonon baths. When an electron is excited to 
the upper Floquet band via a recombination process, it 
quickly “trickles down” via repeated intraband scatter¬ 
ing from acoustic phonons until eventually reaching one 
















FIG. 3. Numerically obtained steady states with radiative 
recombination and coupling to acoustic phonons. Here the 
density is set to half-filling, and we use a 3D acoustic phonon 
bath with huo smaller than the gap A e d ge at the Floquet zone 
edge. The phonon temperature is set to ksT = lCD 2 fif2. We 
keep the phonon and photon densities of states fixed, and only 
vary an overall scale for the coupling matrix elements. The 
full details of the model can be found in Table]!] (a) Distribu¬ 
tion of electrons in the upper Floquet band, F k + = {fl±fk+), 
for several values of k = fcflVV’ rcc / 7r A mter , see Eq. ^ and 
Appendix [C] for definitions. The distributions are fitted to 
a Floquet-Fermi-Dirac distribution at temperature T (solid 
lines). Due to particle-hole symmetry, the distributions of 
holes in the lower Floquet band, 1 — F k _. are identical to 
the distributions shown. Inset Log-Log plot showing the to¬ 
tal density of electrons in the upper Floquet band, n e as a 
function of k. The density n e is normalized to the “thermal 
density” nth = 6.8 x 10“ 4 (see text). The plot demonstrates 
the square root behavior predicted in Eq. Note that for 
large n e , the recombination rates saturate due to Pauli block¬ 
ing. The Floquet band structure is shown in panel (b). 

of the minima of the band. There it joins the Fermi gas 
of excited electrons. Relaxation from the upper to the 
lower band via phonon emission is only substantial near 
the band bottom, where the original valence and conduc¬ 
tion bands are strongly hybridized. The total density of 
excited carriers is determined by a balance between the 
interband excitation and relaxation process. 

As seen in Fig. [3|i, even for relatively large electron- 
phonon coupling strengths the density of excited elec¬ 
trons remains appreciable. As we now explain, this situ¬ 
ation arises from a bottleneck in interband relaxation due 
to the suppression of phonon emission rates for small ex¬ 
citation densities. 

The relaxation bottleneck can be understood by con¬ 
sidering the rate of change of the excitation density 
n e = f fyFfc + of excited electrons in the upper Floquet 
band. In a heuristic model for the regime of low ex¬ 
citation density, recombination transfers electrons from 
the mostly filled states in the “valley” between maxima 
of the lower Floquet band (centered around k = 0), to 


the mostly-empty “hump” in the upper band, providing 
a constant source term for the excitations (see Fig. [l];): 
h T e ec = 7 rec , with 

r k R rjb 

r~J k ( 8 ) 

Here, W[. ec = W k J*~ is the total rate for an electron, 
initially in the lower Floquet band with momentum k, to 
“decay” to the upper Floquet band with any final mo¬ 
mentum (within the constraints of quasi-energy and crys¬ 
tal momentum conservation). Thus is simply the 

recombination rate for a single electron. When we com¬ 
pute 7 rec using W r k c , we take the occupations in the lower 
and upper Floquet bands in the interval —kn<k<kn 
to be one and zero, respectively. It is convenient to de¬ 
fine an average recombination rate in this interval, W , 
whereby Eq. ^ becomes 7 rec ~ (kn/ir)W 1CC . 

Relaxation via interband electron-phonon scattering 
occurs for momenta in narrow regions around ±/cr, from 
the bottom of the upper Floquet band to the top of the 
lower Floquet band. For simplicity, in the discussion 

TTT-k'— TF7^ n ^ er i T77^ n ^ er • 

below we set W k+ = W , where W is an aver¬ 

age value for the transition rates in the active regions 
around ±A:^. The total rate of electrons relaxing from 
the upper to the lower Floquet band is found by sum¬ 
ming the transition rates from occupied states in the 
upper band to empty states in the lower band. The 
corresponding change to the excitation density goes as 
n‘ e nter ^ T J2k' W intC ‘F k+ F k '_. Using particle-hole 
symmetry of the distribution, and Y^k^ k + = Ln e , we 

obtain h I e Ilter « -A inter n^ where A inter = LVF mter . 

The two powers of excitation density appearing in the 
expression for h e come from 1) the density of excited 
electrons available to decay and 2) the density of final 
states available for each electron. When the phonon bath 
is at finite temperature, the picture above is valid when 
n e exceeds the thermal excitation density (see below). 

Importantly, despite the system size L appearing ex¬ 
plicitly in the definition of A lnter , the net relaxation rate 
is in fact system size independent. As explained in Ap¬ 
pendix [bJ the individual rates W kc to scatter between 
specific momentum values k and k' generically scale as 
1/L. The system size independence is restored by the 
increasing number of final states, which scales as L. 

Combining the recombination and interband phonon 
scattering terms, we obtain h e = 7 rec — A lnter rig. The 
condition h e = 0 yields a relation for the steady state 
excitation density n stea dy : 

risteady ~ 1 ^ ^ inter I ■ (9] 

The square root dependence in Eq. (|9| is clearly exhibited 
in our simulation^!, as shown in the inset of Fig. [§. 

Note that in our simulations the bath temperature was 
set to ksT = 0.lAfc H . At this temperature, a “global” 
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Floquet-Fermi-Dirac distribution with its chemical po¬ 
tential set in the middle of the Floquet gap would have a 
small density of excited electrons, n t h (and similarly for 
holes). Here we define the “global” Floquet-Fermi-Dirac 
distribution as a single distribution describing the elec¬ 
tronic occupations in both bands of the system. For very 
low recombination rates, the square root behavior should 
saturate when n e ~ nth- However, throughout the pa¬ 
rameter range used for Fig. UK n e nth, and therefore 
the effect of the finite temperature of the bath on the 
square root behavior is negligible. 

To summarize this section, when radiative recombi¬ 
nation and other Floquet-Umklapp processes are ab¬ 
sent, the system approaches the ideal Floquet insulator 
state (at half filling). Importantly, our analysis shows 
that Floquet-Umklapp processes cannot be ignored: the 
steady state excitation density rises rapidly when the re¬ 
combination rate is increased from zero. In order to fur¬ 
ther reduce the excitation density, additional controls are 
needed. Coupling the system to a Fermi reservoir can 
provide such a control, which we shall study in detail in 
the next section. 


III. COUPLING TO A FERMI RESERVOIR 

In this section we consider the steady state of the sys¬ 
tem upon coupling it to an external fermionic reservoir. 
Our motivation here is twofold: the reservoir serves as an 
additional effective control over the steady state of the 
system, and is a necessary component of transport ex¬ 
periments. However, as we show below, when the driven 
system is coupled to a standard fermionic reservoir with a 
wide bandwidth, photon assisted tunneling significantly 
increases the density of excitations. In addition, even in 
the ballistic regime, phot on assisted tunneling opens ex¬ 
tra channels for transport^ ii | i 2 | 26 | _ w qj discuss how 
such processes can be suppressed using energy filtering, 
thereby allowing for the possibility to control and probe 
the driven system using external fermionic reservoirs. 

In the discussion below we first assume that the dis¬ 
tribution remains homogeneous under coupling to the 
reservoir. This can be approximately satisfied for small 
systems with point-like coupling to a lead, or for sys¬ 
tems where the coupling is extended rather than point¬ 
like. Next we focus on the scenario of a lead coupled at a 
point, where we will discuss the role of inhomogeneities 
and the length scale over which the steady state distri¬ 
bution is controlled by the lead. 


A. Collision integral for a fermionic reservoir 


neling between the reservoir and states of the (undriven) 
system by the Hamiltonian H tnn = Y^t,kv {d' e Cku + 

Cj. „de). The values of the tunneling matrix elements Jt } ka 
depend on the precise forms of the reservoir states {|£)}, 
the Bloch wave functions of the undriven system, and the 
details of the coupling. 

The Floquet states | ipk±) are coupled to the Fermi 
reservoir via the harmonics \(f>k±), as shown in Fig. pb,c. 
These harmonics are spread over a large range of fre¬ 
quencies £k± + nMl. As a result, both the upper and 
lower Floquet bands are coupled to reservoir states in a 
wide range of energies. This coupling is directly mirrored 
in the collision integral for the reservoir. Following the 
spirit of the discussion surrounding Eq. ([6]), we define 
the “bare” rate T^ a for a single electron to tunnel from a 
(filled) reservoir into the Floquet state \ipka{t))-, via the 
harmonic |0£ a ), 

2n 

= X ^ KC«|tftu„K>l 2 s(£ka + nhn - E t ). (10) 

t 


Next, we assume that the reservoir is in equilibrium, with 
the occupation of a state with energy Eg = Ska + nMl 
given by the Fermi-Dirac distribution D{Ef) with chem¬ 
ical potential fx les and temperature T res . To build up the 
integral /J)“ n in the kinetic equation Q, we supplement 
the rates {r£ Q } in Eq. (10) with the occupation factors 
F ka and D(£% a ), with £% a = £ ka + nMl: 


iTa = E [4^(44) - F ka D{£l „)] - (11) 

n 

The first and second terms of Eq. 0 correspond to elec¬ 
trons tunneling into and out of the system, respectively. 


B. Steady state with fermionic and bosonic baths 


How does the coupling to the reservoir influence the 
steady state of the system? The possibility of photon- 
assisted tunneling of particles between the system and 
the reservoir makes the behavior of the driven system 
strikingly different from its equilibrium behavior. 

To understand the effect of the reservoir, it is instruc¬ 
tive to first look at the steady state distribution Fk a of 
Eq. ([4]) in the absence of recombination and electron- 
phonon scattering, I rec = / ph = 0. Staying within the 
homogeneous regime, and setting the left hand side of 
Eq. @ to zero while using Eq. 0 for /^ n , we obtain 


F'kcx. 


E„r^(4J 

E pn 

n ken 


( 12 ) 


The Hamiltonian of the isolated reservoir is given by 
H les = Ee d\dg, where d\ creates an electron in state 
\£) of the reservoir with energy E£. Throughout this work 
we assume that the periodic driving acts only on the sys¬ 
tem, and does not affect the reservoir. We describe tun- 


For a typical metallic reservoir with a wide bandwidth 
(greater than M2), the photon-assisted tunneling rates 
T ka in Eq. ( |Io| may be significant for n ^ 0. Conse¬ 
quently, the sum over n in Eq. (12) leads to steady state 
occupations which differ substantially from those given 
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FIG. 4. Numerically obtained steady states of the system 
coupled to both bosonic and fermionic baths. The top and 
bottom panels show the distributions of electrons in the Flo- 
quet + and - bands, respectively, for increasing strength of the 
coupling to the Fermi reservoir. We characterize the coupling 
strength by the ratio of tunneling and rec ombination rates, 
T = 2r? + /Wr=o ( see Eqs. S and (101 for definitions of 


the rates). Two types of Fermi reservoirs are studied, (a) 
Wide-band Fermi reservoir, whose Fermi level lies in the mid¬ 
dle of the original bandgap (the bandgap of Ha). An increase 
in the coupling strength to such a reservoir leads to a substan¬ 
tial increase in the electron and hole densities n e and rih, due 
to photon assisted tunneling, (b) Energy filtered Fermi reser¬ 
voir, whose Fermi level lies at the resonance energy Eo in the 
original conduction band, i.e., in the middle of the Floquet 
gap of the driven system. The electron and holes densities n e 
and nh are suppressed via the coupling to the narrow-band 
Fermi reservoir. In all panels, the red data points are for a 
half filled system which is disconnected from the Fermi reser¬ 
voir. The other colors correspond to the values of T indicated 
at the bottom. The solid lines are fits to Floquet-Fermi-Dirac 
distributions, with separate chemical potentials for electrons 
and holes in the Floquet + and — bands, respectively. The 
temperature taken for the fits is identical to the phonon and 
reservoir temperature, A?bT = 10 _2 M7. In these simulations, 
the parameters for the photon (recombination) and phonon 
baths were kept fixed at the values yielding the green curve in 
Fig. s while we vary the overall scale of the coupling strength 
to a homogeneously coupled fermionic reservoir. 


by a simple Floquet-Fermi-Dirac distribution Fk a = 

D(£ ka ). 

We now directly illustrate the difficulties which arise 
from coupling the periodically-driven system to a wide¬ 
band fermionic reservoir, in the more general case where 
the system is also subject to electron-phonon coupling 
and radiative recombination, 7 rec , J ph ^ 0. In Fig. [4p, we 
plot steady state distributions for several values of the 
coupling strength to the reservoir. The parameters of 
the bosonic bath (phonons and recombination) are held 
fixed, with values identical to those yielding the green 


(middle) curve of Fig. UK We start at half filling, with 
the chemical potential of the reservoir set in the middle of 
the bandgap of the non-driven system, i.e., we set /z res = 
Eo — MI/2, see Fig. [l] The system-reservoir coupling 
Jt,kv as well as the reservoir density of states are taken 
to be constant 53 . As the reservoir coupling increases, the 
steady state distribution becomes “hotter”, with a higher 
and higher density of excitations. 

The heating effects of the reservoir can be understood 
as follows. In terms of the original (non-driven) bands, 
Fig#, the leading order effect of the reservoir is to pop¬ 
ulate valence-band-like states and to empty conduction- 
band-like states. In terms of the Floquet bands, this 
in particular entails removing electrons from states in 
the lower Floquet band within the momentum window 
—kn < k < kR, and injecting electrons into states of 
the upper Floquet band in the same momentum window. 
Qualitatively, this is similar to the effect of radiative re¬ 
combination, compare to Fig. [3] Strong coupling to the 
reservoir thus leads to a large density of excitations in 
the Floquet bands. To achieve an insulator-like distribu¬ 
tion, as needed for the realization of a Floquet topological 
insulator, these excitations must be suppressed. 


C. Energy filtered fermionic reservoirs 

Interestingly, photon-assisted tunneling can be effec¬ 
tively suppressed if the system-reservoir coupling is medi¬ 
ated through a narrow band of “filter” states (realizations 
are discussed below). For illustration, let us imagine the 
system connected via tunneling to an energy filter: a de¬ 
vice with a finite density of states in a restricted energy 
range, whose states couple strongly to the electron reser¬ 
voir. The filter states hybridize with the reservoir states 
to produce a peak in the continuum density of states 
within the filter energy window. If the system is only cou¬ 
pled to the fermionic reservoir via the filter states, then 
the effective tunneling density of states (TDOS) that en¬ 
ters in the transition rates via Eq. (Tol) quickly falls to 
zero outside the filter window (Fig. [2:). Note that in the 
above discussion we assumed that the filter is not subject 
to the external driv^S 

As a concrete example, consider resonant tunneling 
through a single filter level at energy Filter- Here we 
find tunneling rates with a Lorentzian dependence on en¬ 
ergy: r n ~ 7 / [(£*,« + nHQ - £mter) 2 + Oy/2) 2 ], where 7 
is the level broadening of the filter state due to its cou¬ 
pling to the continuum of reservoir modes. Consider, for 
example, setting Efo ter = I MI = Eq (i.e., at the conduc¬ 
tion band resonance energy). Then, in the limit fit 7 , 
photon-assisted tunneling rates (n 7 ^ 0 ) are strongly sup¬ 
pressed. If the energy filter consists of multiple resonant 
levels connected in series, or a narrow band of states, a 
sharper “box-like” transmission window can be obtained 
(see, e.g., Ref. E3). 

In practice, the energy filter may be realized by cou¬ 
pling the system to the reservoir via a section of large 
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bandgap material hosting a narrow band of impurity 
states deep inside its gap. The intermediate band should 
satisfy three essential requirements: (i) The Fermi level 
should lie inside it, (ii) The band should be narrower than 
the width h£l of a single Floquet zone, as discussed above, 
and (iii) it should be separated from the conduction and 
valence bands of the host material by more than Ml, to 
avoid the direct absorption of photons from the driving 
field. Highly mismatched alloys featuring narrow bands 
of extended states in their bandgaps have been realized 
in the context of intermediate-band solar cells 5 ®^. We 
expect similar methods to allow for the realization of the 
energy filter introduced in this work. Energy filtering 
through quantum dots could provide an alternative ap¬ 
proach. Due to their large size as compared to atoms, 
however, achieving a level spacing exceeding hfl may 
prove challenging (especially at optical frequencies 


D. Steady state with filtered reservoir 


We now investigate how coupling to an energy-filtered 
reservoir affects the steady state of the system. We start 
with the case where phonons and radiative recombination 
are absent, 7 rec = / ph = 0. Throughout the discussion 
below we assume a box-like filter, such that the tunneling 
density of states is strictly zero outside the filter window. 

When the filter window (bandwidth) is narrower than 
MI, photon-assisted processes are suppressed. According 
to Eq. (12), the occupation distribution in the reservoir, 


taken as a function of energy, is then directly mapped 
into the occupation distribution of the driven system, 
taken as a function of quasi-energy (i.e., the occupation 
Fka of each Floquet state | ipka) is determined by a single 
term, D(£ k a ) with fixed n , on the right hand side). In 
the case of the wide-band reservoir, half-filling was en¬ 
sured by placing the chemical potential of the reservoir 
in the middle of the gap of the non-driven system. Here, 
the best choice is to center the filter window around the 
energy of the resonance in either the conduction band 
or the valence band, fi les = ±M2/2, and also to set the 
chemical potential /i res close to the resonance value (for 
the simulations below, we center the filter window around 
the resonance in the conduction band). In this way, the 
reservoir chemical potential will end up inside the Flo¬ 
quet gap. Note that although the chemical potential p, res 
is set to an energy within the conduction the bands of 
the non-driven system, the filtering prevents a large in¬ 
flow or outflow of electrons which would otherwise push 
the density far away from half-filling. Due to the fact 
that this is a highly non-equilibrium situation, however, 
some density shifts away from half-filling are generically 
expected (see discussion below). 

The ideal Floquet insulator distribution can be 
achieved in the situation where the reservoir chemical 
potential is set inside the Floquet gap, and where the fil¬ 
ter window is wide enough to cover the full bandwidth of 
the Floquet-Bloch band structure, but narrower than the 


driving field photon energy hfl such that photon-assisted 
processes are still suppressed. More generically, however, 
the filter window will be narrower than the bandwidth of 
the Floquet-Bloch bands, as depicted in Fig. [TJ-. In this 
case the kinetic equation © with J rec = I ph = 0 does not 
have a unique steady state, as excited electrons and holes 
above and below the filter edges, respectively, have no 
way to relax. However, this is an unstable situation: any 
small scattering rate due to acoustic phonons will allow 
electrons to “trickle down” and fill up all Floquet states 
below the bottom of the filter window. In the absence 
of Floquet-Umklapp processes, the resulting steady state 
will correspond to that of an insulator at finite temper¬ 
ature (assuming the same temperature for the phonons 
and the Fermi reservoir). More specifically, the electronic 
distribution for both bands will be described by a global 
Floquet-Fermi-Dirac distribution with a single chemical 
potential, which is set by that of the reservoir. 

Once we include the contributions of Floquet-Umklapp 
processes such as recombination, the steady state hosts 
densities of excited electrons and holes, n e and n^, re¬ 
spectively (which are generally large compared with the 
thermal density nth)- In the limit of a weakly coupled 
reservoir, the combined density of electrons and holes 
n = n e + nh is determined solely by the recombination 
and phonon scattering rates, as discussed in Sec. |TTj The 
steady state excitation density is further suppressed with 
increasing coupling to the reservoir, as we demonstrate 
below (see Fig. ©>)• 

While the steady state electron and hole excitation 
densities are equal for a half-filled system without cou¬ 
pling to a Fermi reservoir (Sec. 0 . n e and nh need not 
be equal when the reservoir is present, even when the 
chemical potential of the filtered reservoir is placed in 
the middle of the Floquet gap. To see why, note that 
here the Fermi level of the filtered reservoir is aligned 
with the resonance energy in the conduction band 
of the non-driven system. The asymmetric placement of 
the energy window of the reservoir with respect to the 
non-driven band structure generically breaks any effec¬ 
tive particle-hole symmetry, and yields a shift of the to¬ 
tal density away from half-filling. Importantly, the shift 
An = n e — nh can be small, being bounded by n. More 
careful considerations (see Appendix [E]) show that An is 
in fact expected to be significantly lower than n, which 
is confirmed by our numerical simulations (see Fig. [5]d) . 

Staying within the regime of a weakly coupled reser¬ 
voir, let us now consider what happens when the reser¬ 
voir’s Fermi energy is shifted away from the center of the 
Floquet gap. As long as the Fermi level of the reservoir 
remains within the Floquet gap, the occupation factors 
D(£ ka ) in Eq. @ change only weakly, due to the finite 
temperature of the reservoir. Since the rates r£ a are in¬ 
dependent of the occupation of the reservoir, the changes 
in l£“ n are only “thermally activated” by the reservoir’s 
temperature. We therefore expect the steady state of 
the system to be only weakly affected. This implies that 
an interesting situation has been obtained, in which the 
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driven system becomes incompressible , with respect to 
changes of the reservoir’s Fermi level. Once the Fermi 
level enters, say, the upper Floquet band, the density of 
excited electrons in the band is greatly affected. If we 
approximate the distribution of excited electrons by a 
Fermi function, we can expect its Fermi level to track 
the Fermi level of the reservoir. 

As the strength of the coupling to the reservoir is in¬ 
creased, we expect the Fermi reservoir to become more 
dominant in setting the steady state of the system. In 
the limit where the coupling to the reservoir dominates 
all other scattering mechanisms, we expect the steady 
state to be described by a global Floquet-Fermi-Dirac 
distribution, with the same chemical potential as that of 
the filtered reservoir. Note that in this limit, a non-zero 
coupling to the phonon bath is still important in order 
to allow electrons to fill up states from the bottom of the 
lower Floquet band up to the reservoir’s Fermi level. 

The above considerations are confirmed by our numer¬ 
ical simulations, which are given in Figs. [4] and [5j In 
these simulations, we fix the matrix elements describing 
the coupling to the photon (recombination) and acoustic 
phonon baths as in the green (middle) curve of Fig. §b 
and vary the overall scale of the couplings Jg : kv to the 
Fermi reservoir (which are taken to be uniform). The 
reservoir density of states is taken to be constant in 
a window of width Ml /2 placed symmetrically around 
E 0 = I} Ml. The distributions of electrons ( fl a fka ) hr the 
two bands are plotted in Fig. |4j and are separately fitted 
to Floquet-Fermi-Dirac distributions with independent 
chemical potentials fi e and fih for electrons and holes, as 
in Sec. Ill Dl 

Figure [4] clearly demonstrates that for a wide-band 
reservoir, panel (a), the density of excitations increases 
when the coupling to the reservoir is increased; in con¬ 
trast, for a filtered reservoir, panel (b), the density of 
excitations decreases with increasing coupling strength 
to the reservoir. For the filtered reservoir, the chemical 
potential sits at the resonance energy in the conduction 
band, In all fits in Fig. 0 we set the temperatures 

of the Floquet-Fermi-Dirac distributions to be identical 
to the phonon and reservoir temperature. While we ob¬ 
tain excellent fits at weak reservoir coupling, the fits be¬ 
come less accurate when the coupling to the reservoir 
is increased. This arises due to the non-uniform way 
in which the reservoir is coupled to the Floquet bands. 
As in Sec. m we verify that the scattering rates in 
the numerically obtained steady states are significantly 
smaller than the Floquet gap A j. R for all reservoir cou¬ 
pling strengths used (see Appendix [d|) . 

In Fig. [5] we study the densities n e and nh as func¬ 
tions of the strength of the coupling to the reservoir and 
its chemical potential. The figure demonstrates two im¬ 
portant points. First, the steady state densities n e and 
rih are insensitive to small shifts of the Fermi level of 
the reservoir away from the middle of the Floquet gap, 
yielding an “incompressible” behavior dn e ^/dfi rea ~ 0. 
This is demonstrated most clearly by panel (c), which 
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FIG. 5. Electron and holes densities n e and nh in the 
steady state of the system coupled to bosonic baths (acoustic 
phonons and recombination) and an energy filtered fermionic 
reservoir. The figure clearly demonstrates that (1) the steady 
state densities n e and nh are insensitive to small shifts of the 
reservoir’s chemical potential /i res near the middle of the Flo¬ 
quet gap; and (2) a sufficiently strong coupling to the reservoir 
can effectively suppress the electron and hole densities when 
/r r es is within the Floquet gap. Panel (a) shows the total den¬ 
sity n = n e +Uh as a function of the Fermi level of the reservoir 
/ires and the coupling strength ratio T = 2r° + /W£=o- As 
long as /r r es is within the Floquet gap, n e and nh remain low. 
Once /ires enters the Floquet + or — bands, the system be¬ 
comes metallic and the electron (hole) density n e (nh) is set 
by the Fermi level of the reservoir. This behavior is seen in 
panel (b), where we plot n e . To further demonstrate the in¬ 
compressible regime, in (c) we show n as a function of fires for 
several coupling strengths to the reservoir, corresponding to 
the dotted lines in panel (a). Panel (d) gives the the electron 
and hole densities, n e (circles) and nh (squares) for two values 
of fires'- in the middle of the Floquet gap (black) and at the 
edge of the + Floquet band (red). In the first case, the re¬ 
sults explicitly demonstrate the suppression of the excitation 
densities n e and nh with increasing reservoir coupling. Model 
parameters are the same as in Fig. [4] 


shows n vs. fi les (similar plots of A n can be found in 
Appendix [E]) . Second, when the Fermi level of the reser¬ 
voir lies within the Floquet gap, a sufficiently strong cou¬ 
pling to the reservoir can effectively suppress the electron 
and hole densities, giving a steady state close to an ideal 
Floquet insulator. 

The coupling strength at which the reservoir is ex¬ 
pected to significantly affect the steady state excitation 
density can be estimated as follows. Radiative recom¬ 
bination acts to increase the excitation density with the 
rate n T * c = 7 rec defined in Eq. As above, we approx¬ 
imate 7 rec by 7 rec sa (fcj{/7r)W , where kji/7 r represents 
the fraction of states that participate in the recombina¬ 
tion process and W is the average recombination rate 
in this interval. Meanwhile, the reservoir can extract ex- 
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citations at a rate h * un ss —T°n e , where r° = + 

characterizes the rate for a single electron to tunnel in 
or out of the system. When the reservoir is weakly cou¬ 
pled, the steady state excitation density is controlled by 
the rate of phonon-mediated interband relaxation, as dis¬ 
cussed around Eq. ([8]). The reservoir begins to play an 
important role when T 0 n e (/c/j/7r)W , or equivalently 

when 

TE(2r°/wa>WK). (is) 

This relation is indeed borne out in Fig. [5jl. where 
n e ss 0.05 and kf i/tt « 0.3, and the excitation density is 
suppressed for log 10 T > 1. Note that Fig. 5jl also shows 
a small non-zero value of An <C n when the Fermi level 
is in the middle of the Floquet gap (black symbols in 
Fig#), arising from the asymmetry of the coupling of 
the reservoir to the two Floquet bands. 


E. Point coupling to a lead 


In many experimentally relevant situations, the system 
is coupled to a lead at a single point. What is the spatial 
dependence of the steady state in this situation? So far 
we have discussed homogeneous steady state distribution 
functions {Fk a }, which depend on momentum and band 
indices but not on position. A homogeneous steady state 
can arise in several situations. For the bosonic baths dis¬ 
cussed earlier, we assumed a uniform coupling through¬ 
out the system. Therefore by themselves the bosonic 
baths yield a spatially homogeneous steady state distri¬ 
bution. Clearly, if in addition we introduce a fermionic 
reservoir which is coupled homogeneously throughout the 
system, a spatially homogeneous steady state is expected. 
In addition, for a lead coupled at a single point, there are 
still two limits in which the steady state remains uniform: 
(1) absence of any other sources of dissipation; and (2) 
the limit of a small system size. In the latter case, a uni¬ 
form distribution is obtained when the level spacing of 
the system’s single particle states is larger than the tun¬ 
neling rate to and from the reservoir; this corresponds to 
a tunneling time that exceeds the time required for an 
electron to traverse the system. 

For larger system sizes, where the stringent criterion 
above is not met, the steady state need not be spatially 
homogeneous. If the tunneling rates are comparable to or 
larger than the level spacing, the coupling to the reser¬ 
voir can yield nonzero values of the “coherence” terms 
(fl a fk'/ 3 ), which generically cause spatial inhomogeneity. 
Calculating the full set of such coherences is a formidable 
task. Fortunately, we can gain an intuitive understand¬ 
ing of the form of the inhomogeneous steady state by 
considering the dynamics of a spatially-dependent exci¬ 
tation densities n e (x,t) and nh(x,t). Close to the lead, 
placed at x = 0, the excitation density will be affected 
by the lead and will roughly correspond to that found for 
a homogeneous system-reservoir coupling. Far from the 


lead, we expect the excitation densities to relax to bulk 
values 7ig Ulk and n]? ulk . Below we estimate the “healing 
length” over which this transition occurs. 

Due to fast intraband electron-phonon scattering 
(which is still slow compared with the driving frequency 
and the on-resonance Rabi frequency A k R /h), carrier mo¬ 
tion on time scales much larger than the driving period 
is expected to be diffusive. The corresponding diffusion 
constant can be estimated as D = v 2 /W mtra ', where v 
is a typical velocity of the excitations 59 and W nntra is a 
typical intraband scattering rate from acoustic phonons, 
both taken in the steady state. Focusing on the situ¬ 
ation near half-filling and incorporating the source and 
sink terms due to recombination and interband phonon 
scattering discussed in Sec. [TTJ we obtain two reaction- 
diffusion equations for the electron and hole densities, 
d t n\ = D d 2 n\+Y ec - A intel n e n h , with A = e,h. Adding 
and subtracting these equations, we find the reaction- 
diffusion equations governing the total and offset densi¬ 
ties n and An, 

d t n = Dd 2 x h + 2 7 rec - iA inter (n 2 - An 2 ), 
d t An = D d 2 An. (14) 


For the boundary conditions for the above equations, 
we use A n(x = 0) and h{x = 0) which are set by the 
lead, as well as d x An = d x h = 0 for x 0, which 
corresponds to no net flow of particles into the system. 
Eq. (14) entails two main consequences for the spatial 
distribution of the steady state, d t n a (x,t) = 0. First, 
the shift of the total density of electrons from half filling, 
An, is in fact homogeneous across the system, and set by 
the lead. Furthermore, linearizing Eq. (14) around the 
bulk steady state gives the healing length 


/ Dn bulk \ 1/2 
\ 47 rec ) ’ 


(15) 


where we can approximate nbuik = 2y £ 7 ?>V 1CC / 7r A mter 
from Eq. ©• Here we neglect corrections due to a small 
An. For system sizes smaller than £, a lead coupled at a 
point can be effective in setting the distribution through¬ 
out the system. In such a system, for a sufficiently strong 
point coupling to a filtered lead, a Floquet insulator dis¬ 
tribution can be achieved, as was shown for the homoge¬ 
neous case in Sec. IIII Dl 


IV. SUMMARY AND DISCUSSION 

The ability to control and probe non equilibrium quan¬ 
tum many body systems poses one of the most outstand¬ 
ing challenges in modern condensed matter physics. In 
this paper, we analyzed steady states in a model for a pe¬ 
riodically driven semiconductor, and demonstrated the 
means through which these steady states can be con¬ 
trolled. We considered the open system dynamics of a 
resonantly driven electronic system coupled to acoustic 
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phonons and the electromagnetic environment, as well as 
an external fermionic reservoir. The couplings to these 
baths have two complementary roles: they allow energy 
relaxation, but may also induce processes which lead to 
heating. Motivated by the prospect of realizing Floquet 
topological insulators, our goal was to find the conditions 
under which the steady state resembles a band insula¬ 
tor. Importantly, we focused on the regime where the 
scattering rates in the steady state are smaller than the 
drive-induced Floquet gap. Only in this regime could we 
expect to observe effects requiring quantum coherence, 
such as drive-induced topological phenomena. 

Starting with the case where the system is coupled only 
to the bosonic baths, we have shown that the system can 
approach a Floquet insulator steady state with an added 
density of excitations in the two Floquet bands. The 
density of excitations is controlled by the ratio of radia¬ 
tive recombination and electron-phonon scattering rates, 
and can be small for experimentally-relevant parameter 
values. We found a square root dependence of the excita¬ 
tion density on the above ratio, see Eq. ([9 ), which implies 
that additional controls are needed to fully suppress the 
deviations from the Floquet insulator state. Next, we 
considered the effects of coupling to an external Fermi 
reservoir, which plays two important roles in our setup. 
First, the reservoir is a crucial component for transport 
experiments. Importantly, we show this coupling signif¬ 
icantly increases the density of excitations, unless the 
reservoir is coupled through an energy filter. Second, 
the energy-filtered reservoir can serve as an additional 
control to reduce the density of excitations, bringing the 
system closer to the ideal Floquet insulator state. 

A main motivation for our work is the prospect of ob¬ 
taining a Floquet topological insulator: a driven system 
with an insulating bulk but with conductive edge and 
surface modes. What are the implications of our results 
for transport? Interestingly, our results show that even in 
the presence of a finite density of excitations, the steady 
state of the driven system can exhibit an “incompress¬ 
ible” behavior: the steady state is unaffected by small 
changes in the chemical potential of the energy filtered 
reservoir, as long as it situated near the middle of the 
Floquet gap®. This behavior is shown in Fig. [ 5 J In ad¬ 
dition, when energy filtered leads are used for transport, 
photon-assisted conduction channels are suppressed. In 
the case of neutral particles, the incompressibility im¬ 
plies insulating behavior: no current would result from a 
small source-drain bias between two spatially-separated 
energy-filtered leads. This follows from the insensitivity 
of the steady state to the leads’ chemical potentials. 

When we consider charged particles, however, the non¬ 
zero density of excited carriers is expected to yield a fi¬ 
nite resistivity even in the paramter range correspond¬ 
ing to the insulating regime above. Given the diffu¬ 
sive nature of particle motion in the system, we esti¬ 
mate the local resistivity a^ x (x) using the Drude form: 
cr~ 1 (x) = |m*|l / F mtra /[e 2 n(a;)], where m* is the effec¬ 
tive mass around the Floquet upper (lower) band mini¬ 


mum (maximum), e is the electric charge. Consider now 
a two-terminal transport measurement using energy fil¬ 
tered leads, through such a system of charged carriers. 
If the system is small enough such that the steady state 
is spatially homogeneous, a sufficiently strong coupling 
to energy filtered leads can, in fact, suppress the den¬ 
sity of excited carriers and yield nearly-insulating be¬ 
havior. For larger systems, with a spatially inhomo¬ 
geneous steady state, the total resistance R is given 
by the sum of series resistances, R = / Q L dxa~ 1 (x). 
The bulk of the system gives an extensive contribution 
-Rbuik ~ L\m.*,\W lntla / [e 2 n bulk ]. Interestingly, if the lead 
coupling is strong, the excitation density near the ends 
of the system will become very small and thus give a 
large contribution R en( j to the resistance. Therefore, for 
a fixed system size, the system may obtain insulating be¬ 
havior in a two-terminal measurement upon increasing 
the coupling strength of the lead. 

The analysis of steady states in driven electronic sys¬ 
tems is currently the subject of intense activity (see e.g., 
Refs. 1251 and Hill . In Ref. SS) time evolution and steady 
states after a quench were studied for a 2D semiconduc¬ 
tor with a topological Floquet spectrum. There, the au¬ 
thors considered open system dynamics with momentum 
conserving interactions with a bosonic bath, and found 
regimes exhibiting quantized Hall conductivity. In addi¬ 
tion, Ref. ED studied a resonantly-driven electronic sys¬ 
tem where the only relaxation pathway was through an 
external fermionic lead, and found that a grand canonical 
distribution could be obtained under finely-tuned condi¬ 
tions. In our work, which was carried out in parallel, we 
considered the combined effects of momentum and en¬ 
ergy relaxation through the coupling to acoustic phonons 
and the coupling to an external fermionic reservoir. No¬ 
tably, we included the inevitable effects of heating due 
to Floquet-Umklapp processes, as exemplified by radia¬ 
tive recombination. Importantly, momentum relaxation 
played a crucial role in establishing the Floquet insulator 
steady state under these conditions. 

To make a connection with experimentally relevant 
regimes, we compare our model parameters with those 
accessible in solid state systems. Consider a drive fre¬ 
quency of = 2n x 100 THz, which translates to 0.4 
eV in energy units. Correspondingly, the parameters 
used in Sec. [TTJ yield a Rabi frequency (Floquet gap) 
of Ak R /h ~ 27t x 10 THz (translating to 40 meV), 
and a characteristic phonon relaxation time scale of 
Efc' W* + + ] _1 = 1 ps. This is the total relaxation rate 
out of the state k = 0 in the upper band. The recombi¬ 
nation time scales, (Wj^n)” 1 used to obtain the steady 
state distribution in Fig. [3j are [1 /zs, 60 ns, 3 ns, 180 
ps, 10 ps]. When coupling to the fermionic reservoir was 
introduced in Sec. m we fixed the recombination time 
scale at 3 ns; the steady state distributions in Fig. [4] cor¬ 
respond to tunneling times (Ffc = o,-)^ 1 of approximately 
[200 ps, 30 ps, 3 ps]. In Fig. [5| the tunneling time from 
the reservoir varies from 30 ns to 3 ps. Note that these 
values are in line with those in typical semiconductor 



15 


nanostructures, where tunneling times can vary widely 50 . 

Although our model is inspired by resonantly-driven 
semiconductors, we expect our conclusions and formal¬ 
ism to be relevant to a broad variety of driven-dissipative 
systems including cold atomic gases. Our results also 
have important implications for Floquet topological in¬ 
sulators. Indeed, they provide a roadmap towards the 
practical realization of the Floquet insulator state, which 
is key to observing quantized transport in Floquet topo¬ 
logical insulators. We expect engineered reservoirs of 
carriers to be particularly useful in this context, allow¬ 
ing to perform transport measurements while stabilizing 
insulating-like steady states. 

Several aspects of the problem require further study. 
In this work we have not addressed the effect of inter 
particle interactions. Floquet-Umklapp processes involv¬ 
ing inter-particle scattering give an additional channel for 
the system to absorb energy from the driving field and 
increase the number of excitations. However, our current 
work demonstrates that coupling to a bath of phonons, 
can keep the heating and excitation density under con¬ 
trol. Another important direction is a careful study of 
the inhomogeneous steady states of Floquet topological 
insulators, which is crucial in order to predict the edge 
and surface response of these systems. 
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Appendix A: Kinetic equation in the Floquet basis 

In this appendix we discuss the key points in the 
derivation of the Floquet kinetic equation, represented 
schematically in Eq. Q of the main text. Specifically, 
the aim is to derive a system of differential equations 
which describe the time evolution of the Floquet state 
occupation factors F ka (t) = {ft a (t)f ka (t)), where f ka (t) 
and f ka (t) are the creation and annihilation operators for 
Floquet states defined above Eq. The time derivative 
Fka = couples to an infinite hierarchy of 

higher and higher order correlation functions. The main 
approximation is to truncate this hierarchy at the lowest 
non-trivial order, and obtain a closed system of evolution 
equations. We now outline the required steps. 

1. Basis transformation and dressed matrix 
elements 

The first important step is to express the electronic 
terms in the Hamiltonian in terms of the Floquet creation 
and annihilation operators. Using Eq. |3|, the transfor¬ 
mation is made via 

4 = EE e<(£fc “ +nn)t ^ai^>/Lw. 

a n 

c k „ = E E e-^ +mn)t (vk\Wp)fk 0 (t), (Al) 

P m 

where \vk) is the Bloch function corresponding to crystal 
momentum k in band v of the non-driven system, and 
f ka (t) and f k a (t) are creation and annihilation opera¬ 
tors for the Floquet state \if ka {t)) (we simplify notations 
by setting H = 1 here and everywhere below). Using 
these relations, we write the electron-boson interaction 
and system-reservoir tunneling Hamiltonians as 


Hint = E E E q x )ft_ g ^ a ,(t)f ka (t)(b_ q + 6+), 

kq aa.' n 






k£ n 

with “dressed” matrix elements 

«*) = E E k - q*)Wk\W a ), 

v m 

J$a = Y, J eM<t>: k Wk). 


(A2) 

(A3) 


(A4) 

(A5) 


Here we use Q and J to indicate the coupling matrix 
elements in the Floquet basis. 


Note that in Eq. (A4) we have imposed lattice mo¬ 
mentum conservation of the electron-phonon interaction, 
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Gi,f (q x ) oc Sq^^k-k 1 - Equations (|A2|) and (|A4|) arise from 


Eq. 0 of the main text. Sets of equivalent scattering 
processes can be identified based on the following useful 
relation: 


q x ) = [G { ao ?\k + q x , -q x )]*. (A6) 

The explicit form of the dressed matrix elements shows 
that the overlaps {<t>™ a \kv) between the original states 
and the Fourier components \(/>™ a ) of the Floquet modes 
are crucial in determining the rates of the different Flo¬ 
quet scattering processes, as discussed in the main text. 


/L(f), Similar relations hold for the annihilation op¬ 
erators. Let U(t,t') be the single particle time evolu¬ 
tion operator corresponding to the Schrodinger equa¬ 
tion ;), with the Hamiltonian H(t) = 

Y^k c lv H ™'(k,t)c kl/ >, where H(k,t ) = H 0 (k) + V(t) is 
defined in Eq. (Ilj) in the text. The operator f ka (t) satis- 
fies fl a (t) = U(t, t') f ka (t')W (t , t'), which can be written 
in the differential form: 


= {H(t),fl a (t)\. (A7) 


Matrix elements for Floquet- Umklapp processes 

An interesting situation occurs when the coupling to 
the driving field is defined by a vector g [see Eq. 0] such 
that <j|| = 0 [see Eq. ([2j) for the definition of g]. This com¬ 
monly occurs in experimentally relevant materials driven 
by optical fields. Here, the Fourier harmonics | (j)™ a ) have 
a fixed band character for m of fixed parity: for example, 
in the convention used throughout the paper and set be¬ 
low Eq. (3), \(j)™ a ) is proportional to | kv) for m odd and 
proportional to \kc) for m even, see Fig. 2a. Additionally, 
note that scattering by phonons preserves the band char¬ 
acter, G uv ’ oc . As a consequence, under these condi¬ 
tions Floquet-Umklapp processes involving phonons are 
forbidden for n odd. Therefore, if in addition to g» = 0 
the phonon bandwidth is less than the driving frequency 
Cl, all Floquet-Umklapp processes, including both even 
and odd n, are not allowed. 


2. Equations of motion 

We now study the equations of motion for the Flo¬ 
quet state populations F ka = ( fl a (t)fka(t ))■ The popu¬ 
lations are the diagonal part of the “polarization matrix” 
p£?{t) = (/fc'a'( i )/fcc(t))- In addition to the popula¬ 
tions, this matrix also characterizes coherence between 
Floquet states with different crystal momenta and/or 
band indices. This off-diagonal part may be important 
for the dynamics and for characterizing steady states. In 
the main text we focus on steady states in a regime where 
the off-diagonal part of the polarization matrix can be ne¬ 
glected. Here we derive the kinetic equation in a more 
general context, including the full polarization matrix, 
and discuss when and how the off-diagonal parts may be 
neglected. 

As a preliminary, we note the following impor¬ 
tant property of the Floquet state creation operators 


This expression will be used below. 

The derivation of the kinetic equation proceeds along 
standard lines, as explained in detail in, e.g., Ref.Hfr, The 
main difference from the usual case (i.e., for non-driven 
systems) is the appearance of the “dressed” matrix ele¬ 
ments in the interaction Hamiltonians. Below we set up 
the calculation and point out where these terms appear, 
and where special considerations are needed to complete 
the derivation for the case of a periodically driven system. 

We seek the time evolution of the Floquet state pop¬ 
ulations Fk a (t). However, since these populations are 
special cases of the polarizations P(t) defined above, 
for k = k ', a = a ', we begin with the more general ex¬ 
pression for the time derivative of (t): 

idt(fl a ,mka(t)) = (A8) 


where H{t) is the full single particle Hamiltonian (includ¬ 
ing driving) as defined above, and H tot = H(t) + F [b + 
Hint + H res + H tun is the total Hamiltonian including the 
baths and the system-bath coupling. The commutator 
in Eq. (A8) includes two types of contributions, arising 
from: (1) the time derivative acting on the state with 
respect to which the average is taken, and (2) from the 
explicit time dependence of the op erat ors f' k i a i(t) and 
fka(t). The latter are given by Eq. (A7) and its Hermi- 
tian conjugate. 


To simplify the expressions below, we introduce a more 
compact notation in which k and the Floquet band in¬ 
dex a are compressed into a single index a. In this no¬ 
tation, the dressed electron-phonon coupling matrix el¬ 
ements will be wri tten as Q^ a {k,q x ) = G^l(q x )- The 
commutator in Eq. ( |A8[ ) has two non-trivial terms related 
to the system-boson and system-reservoir couplings Ft -, n t 
and Htun, [fl(t)fb(t), H int \ and H tun ], respec¬ 

tively. The system-boson coupling produces the following 
contribution: 
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mt)m,H int \) = j2E ei(£b ~ £a ' )teinnt ^(^ (fi(t)fa'(t)(b„ q + bD) 

a' q n 

~ ^^ e ^'-^ )t e“^i(fe)(/i,(t)/ 6 (t)(6- q + b\)), (A9) 

a' q n 


while the system-reservoir coupling leads to 

([fl(t)f b (t),H ies ]) =Y.Y. ei£bteinatj $ (faWe) 


t n 


~T,T, e ~ i£ate ~ inntj ta* (fb(t)d\). 


(A10) 


t n 


Note the appearance of “mixed” correlators such as 
{fl{t)fa'(t)(b- q + b\)) and ( fl(t)d e ) involving both sys¬ 
tem and bath degrees of freedom, which appear in 
Eqs. (A9) and ( A10| . The expressions are very similar 
to those that would be obtained for a non-driven system, 
except that here we find an additional sum over n which 
accounts for the harmonic structure of the Floquet state 
wave functions. 

In order to describe scattering between Floquet states, 
we need to solve for the equations of motion of these 
three-point correlators. To do so, we must evaluate ex¬ 
pressions such as 

id t {fi(t)f b {t)b_ q ) = ([ fl(t)f b {t)b_ q ,H tot - H(t)]} 
idt(fl(t)de) = - H(t)}). (All) 

Similar expressions are also needed for fb(t)b\) 

and id t (fb(t)d\). 


and ( fifbd\d e >). 


The commutators in Eq. (All) generate many terms. 
The corresponding calculation is straightforward, but 
somewhat tedious. As above, the primarydifference from 
the textbook case of a non-driven system^ is the appear¬ 
ance of sums over Floquet harmonic indices. 

Mathematically, the crucial point is that the commu¬ 
tators in Eq. (All) give rise to higher order correla- 


I 

tion functions such as {flfthfdi>\b q ,, \ Ja Jb^ 

In the first case we split the averages into prod¬ 
ucts of averages of fermionic and bosonic bilinear 
operators: (flflfbfdb\b q >) « ( flfd)(flfb)(b\b q ') - 
(flfb)(flfd)(b\b ql ), etc. The fermionic averages involv¬ 
ing system operators just give the polarizations P“ de¬ 
fined above. We take the averages of the bosonic op¬ 
erators with respect to a thermal distribution with in¬ 
verse temperature /?: (b^ q b q :) = 5 qq 'N(uj q ), where N(e) = 
1/(1 — e -/3E ) and uj q is the frequency of bosonic mode q. 
Likewise, we split the averages involving reservoir degrees 
of freedom as (flf b d],de') ~ (flf b )(d\de')- For the Fermi 
reservoir, we take (d\d^) = See'D(Ee), where D(E ) is the 
Fermi-Dirac function with temperature T and chemical 
potential q les . For brevity, below we use De = D(Ee). 

Through the above approximations we close the equa¬ 
tion of motion hierarchy. After splitting the averages 
on the right hand sides of Eq. (All), we integrate 


them from time 0 to t to find the correlation functions 
( fa(t)fb(.t)b q }(t ) and (fl(t)de)(t) needed as input for the 
equations of motion of the polarizations, Eqs. (A8), (A9) 
and ( |A10[ ). 

To give an explicit example, we focus on one term 
which arises from the system-reservoir coupling, 


idt(flde) = E e (p a de) + J2J2 e~ l£bt e~ inQt J$* [P b a ( 1 - D t ) - (, 5 ab - P b a )D e } . (A12) 

b n 

The calculation for other terms yields similar expressions. A straightforward formal integration, taking (f^d b )(t = 
0) = 0, then yields 


(fide) = -VV [ dze-^ £b - E ^ e - m ° s [P b a (l - D e ) - ( S ab - P^)D e ] . (A13) 

1 „ J o 


b n 


The next step is to introduce this result and its coun¬ 


terpart for ( d\f a ) into Eq. (A10), for the contribution 
of the system-reserv oir c oupling to the evolution of the 
population P“ , Eq. (A8). Doing so, we obtain: 
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dflfa', JTtun]} = - E E f dze -^-Zt)* e -imn Z [pa (1 _ ^ _ pa )£)/ ] 

* a mu J o 

--VV j{n)*j{rn) e -i { S a - Ee )t e ~inQ.t f dze i{£ b -E e )z jmQz _ pbjp^ _ pb^ _ ( Al4 ) 

1 ou ’ " -'c 


£6 ran 


r 


Importantly, notice that the right hand side of 
Eq. (A14) couples the evolution of the population F a = 


P“ to both the diagonal and off-diagonal polarizations 
P h “. Thus in principle we do not have a closed set of 
equations for the populations alone. In particular, for 
transient behavior (e.g., at early times when the driving 
is just switched on) such terms can not be ignored. 

Close to the steady state, we may expect the off- 
diagonal polarizations (coherences) to be small under cer¬ 
tain circumstances. For a homogeneous system where 
the steady state maintains translational invariance, the 
polarizations in the steady state are diagonal in the elec¬ 
tronic crystal momentum, P kc ^ oc 5 k k'■ Furthermore, 
coherences between the two Floquet bands can be sup¬ 
pressed in the steady state under suitable conditions, 
which are discussed at length in Appendix [D] These con¬ 
ditions are expected to be met for weak system-bath cou¬ 
pling, and we have verified that the steady states result¬ 
ing from our simulations are indeed in this regime (see 
Appendix [D]). 

For strong system-bath coupling, the conditions dis¬ 
cussed in Appendix [D] might not be met, and a more 
complicated situation may arise. There, the particular 


form of system-bath coupling may try to drive the sys¬ 
tem towards specific states other than the Floquet states. 
For example, relaxation may occur into the eigenstates of 
the non-driven system. The competition between driving 
and relaxation may then lead to steady states featuring 
significant inter-Floquet-band coherences. 

In this work we focus on the case of homogeneous 
steady states, with weak (but nonetheless realistic) sys¬ 
tem bath coupling. We negl ect al l off-diagonal coher¬ 
ences, setting a <5 a & in Eq. (A14) and similarly for all 
other terms in the equations of motion. Additionally, in 
the sums over Fourier harmonics we only keep the terms 
for which n = m; when the Floquet state populations 
evolve slowly on the timescale of the driving period, the 
terms with n ^ m give rise to fast oscillations and thus 
produce negligible contributions. With these two simpli¬ 
fications, the standard Markovian approximation yields 
the full Floquet kinetic equation: 


9tF ka = It 


rrec 

1 kot 


L ka ? 


(A15) 


with the collision integral for electron-phonon scattering 
given by 


t = ?EE q x ot' Fka,N {?VjJ q) FfcaFfc—q x ot' (1 T S(^Skct Sk—q x ot' fudq TlfrFV) (A16) 


a.' q n 


+ ?EE \fil q x a.'Fkoi (1 N{fVjJq)') FkaFk—q x a'N(hjjJq)^ 8{Ek—q x ot' &kot ftjjJq ~"b Tlfa&V). 


a'q n 


Tunneling in and out of the Fermi reservoir is described by (see Eqs. (|10[) and (11) in the main text): 


27r 


^“ = tEE \J$L\ I * 3 [FkaD{E t ) - F ka ( 1 - D(E e ))\ S(£ ka -E e + nfifi). 


H n 


(A17) 


The collision integral corresponding to radiative recom¬ 
bination looks identical to that for electron-phonon scat¬ 


tering in Eq. (A16), with the matrix elements q x ) 


replaced by the appropriate ones for coupling to the elec¬ 
tromagnetic environment. In our model, the matrix ele¬ 
ment for coupling to bath photons is purely off-diagonal 
in the basis of the conduction and valence bands of the 
non-driven system. This model is motivated by the form 
of radiative transitions for electrons near k = 0 in many 
experimentally relevant materials. For simplicity, we 
modeled recombination as “vertical” transitions^, giv¬ 
ing G rec = 5 rec (l - <W)<S?*,o4,k'- 


I 

According to our convention in Eq. ([3]), the fact that 
G rec oc (1 — 4„/) requires a coupling between \<j>™ a ) and 
where m an d m' are separated by an odd inte¬ 
ger for the case g» = 0. Furthermore, the conservation 


of energy expressed by the delta function in Eq. (A16) 
requires n to be negative. Therefore, in our model re¬ 


combination only acts through terms in Eq. (A16) with 


n < 0 odd. The dominant contribution comes for n = — 1 
for weak driving. Correspondingly, the emitted photon 
energy is large (on the order of the driving frequency), 
and hence we set all Bose occupation factors for photons 
to zero (i.e., only spontaneous emission is included). 
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Finally, to get the collision integral (A161 into the form 
of Eq. 0 in the text, we integrate over the delta function 


to get the density of states for bosons with momentum 
q x parallel to the system. This gives 


£ = f££i#l(* ,<a)| 2 [Fk'a'F ka N(A£ n ) - F ka F k ’ a f(l + N(A£ n ))] p qx ( A£ n ) (A18) 

k'a' n 

+ f££ i#>(* ,Qx)\ 2 [F k , a ,F ka (l + N(-AEn)) - F ka F k i a ’N(— A£„)] p qa .(—A£ n ), 

k'a' n 


where in the above q x = k — k' and A£ n = £ ka — £ k ' a ' — nhFl. 


Appendix B: System size scaling of transition rates 


In this section we discuss the scaling of the electron- 
boson scattering rates W £ a “ with the system size. We 
will show that in the limit of a large system, the rates 
scale as ~ 1/L. As we explain below, this implies that 

both 7 rec and A inter = LW' nte \ defined in Eq. ^ and the 
discussion below, are independent of system size. An im¬ 
portant consequence is that the excitation density n ste ady 
is also independent of the system size, as one would nat¬ 
urally expect. We first focus our discussion on the radia¬ 
tive recombination, i.e. interaction with a photon bath, 
and then explain how it can be easily applied also to a 
bath of phonons. To simplify the discussion, we illustrate 
the scaling using a non-driven toy model, but the discus¬ 
sion can be easily generalized for transition rate between 
Floquet states in a driven system. 

We consider an electronic Bloch Hamiltonian of the 
form H 0 (k) = [2A(\ — co$,{ka)) + E g&p ]a z . We define 
the Bloch states as | ka) = Si=o ^ elkx \x , a), where 

c, v correspond to the positive and negative eigenvalues of 
a z , a is the lattice constant and L = Na is the electronic 
system size. 

The electron photon interaction Hamiltonian, in the 
rotating wave approximation, is given by F [ lnt = 
T,q H mt{q), with 


x 

(Bi) 

where in the above cj. a are creation and annihilation op¬ 
erators for Wannier states in the conduction and valence 
bands, and \M q \ depends on the volume of the electro¬ 
magnetic environment as \M q \ ~ l/y/V env . Note that 
Hi n t(q) is diagonal in the lattice coordinate x. The rate 
for recombination from \k,c) to | k',v) is then given by 


Wj? c v 


27r 


£ 


^£- 


,i(k—k' -\-q)x 

N 


2 

3(F kc E k ' v hidq ), 


(B2) 

where { E ka } are the eigenenergies of F[ 0 (k). Impor¬ 
tantly, the photon momentum lives on a different recip¬ 
rocal lattice than the momenta of the electronic system, 
q = f^-n. For simplicity, we drop the q dependence 


of M q . Summing over the transverse photon momenta 
q l yields a 2D density of states for the transverse modes 
with q x held fixed, 




\M\ 2 

n 

q x 


£ 


a i(k—k'+q)x 


Pq x (Eke Ek'v')* 


(B3) 

Note that p qx (E kc — E k i v ) has dimensions of En * and 
scales as L 2 nv . The photons emitted by the radia¬ 
tive recombination transition have a typical energy of 
E gapi and therefore the corresponding photon momen¬ 
tum Tig* = F/gap /c plays an important role in the calcu¬ 
lation of the rates. For simplicity, we set the density of 
states p qx {E kc —E k > v ) to be a constant p^ D ^ for \q x \ < g*, 
and zero otherwise. This gives 


w£? 


2tt M 2 p {2D) 
H N 2 


£ 


q x =-q* 


2 


"y ^ — k'+q)x 

x 


(B4) 


Assuming a large environment volume, we can write 


k , v= 2_vM'X 

" kc t- 


2 J? D )r 


N 2 


f 1 - dq 

1 - e iqL 

l-q, 27T 

1 _ e i(k — k'+q)a 

( 


Recalling that M ~ 1 / y/V env we see that the factor 


M 2 p {2D) L R 


is independent of environment size. 


The calculation now amounts to evaluating the inte¬ 


gral in Eq. (B5). Using integers to represent momenta 
as in k = qfn, we denote this integral by g^v(n — n'), 
where the N subscript denotes the fact that the integral 
depends on the system size L = Na. We are interested 
in the scaling of this integral with N. We define the di¬ 
mensionless variable q = qL 1 and divide by TV 2 for later 
convenience, whereby the integral becomes 


fliv(m) 

N 2 


1 

L I 


q * Na d ~ sm 2 (q/2) 

q,Na q N 2 sin 2 ( 2 ^ [2ttto + g]) 


(B6) 


Note that in the prefactor on the right hand side above, 
L is the electronic system size. In order for W k c v ~ 1/L, 
which guarantees that e.g. ro s teady remains independent 
of system size, one has to have that ~ 1/L. 

In the following, we assume g*a = Ea ff a -C 1, which 
means the photon wavelength is much larger then the 
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lattice spacing of the system. We furthermore consider 
the limit where the system size is larger than the photon 
wavelength, N l/(<?*a). We start again from Eq. (B 6 ). 
Clearly, for gjq{m)/N 2 to be of order 1/L, we must have 
27 t|to| < Nq*a, which guarantees that the integral picks 
the contribution where the sin function in the denomina¬ 


tor of (B 6 ) vanishes. Physically, this corresponds to the 


requirement that | k — k' | < g*. 

We now need to check how the integral in Eq. (B 6 ) 
scales with N. We do this explicitly for m = 
0; the result can be generalized for any 27 t|to| < 

Nq*a. We break the integral into three integration re¬ 
gions: (1) [-^q*Na, y/q*Na\ (2) [y/q*Na, q*Na] and 
(3) [— q*Na, —y/q*Na\. In region (2), we can give an 
upper bound to the integral by 


nq*Na 


dq 


N—>oo 


N 2 sin 2 ( \jq*aj4N) 


>C, 


(B7) 


where C is a constant. The same result applies to the 
integral in region (3). In region (1), we expand the de¬ 
nominator to obtain 


fV^N-a sin 2 (g72) 

/ dq -— 

J — \/q*N a 4 Q 


1 ---h 

12iV 2 


(B 8 ) 


The first term in the above expansion clearly gives an 
order 1 contribution, while the rest of the terms vanish 
in the limit of large N. 

The result of the above analysis is that < 7 at (to) ~ 
g(m ), where g{m) is independent of system size. The 
full expression for gi\r(m ) may contain terms that scale 
slower than N 2 /L with the system size. Finally, insert¬ 
ing this result this back into Eq. (B5), we arrive at the 


scaling W kv c ~ 1/L of rates with the system size, in the 
limit of a large system. 

Putting this into the definition of the total rate of re¬ 
combination out of the state | fc, c), defined by W/ ec = 
Zk>W& v , we get 


wr c = L 


dk w 


27r 


kv 

kc 


E, 


9 a P \\fkv 

nhc kc 


(B9) 


Therefore WjC ec is independent of system size, as 
promised. Likewise, the total rate density, j Tec = 
T zCfc fc' Wfc = f || W k ec is independent of system size. 

Note that the factor in Eq. (B9) accounts for the 


photon density of states in the longitudinal direction, 
whereby p 3D = q^?fp 2D - 


k' + q = k. We start from the analogue of Eq. (B2), for 
phonon scattering rates. Performing the sum over x we 
obtain a factor of N 2 S k ’_ kqx - In this case, however, the 
factor M 2 Pq D scales as ~ 1/L. Therefore, W k ,° ~ 1/L. 


Appendix C: Numerical simulations 

In our numerical simulations, the steady state distri¬ 
butions were obtained by direct evolution of Eq. @- 
The results are independent of the initial distribution 
{Fkait = 0)}. 

In Sec. [TT] we discussed the square root depen¬ 
dence of the excitation density n ste ady on the ratio 
7rW ieC /fcRlE lnter . This behaviour was observed in our 
numerical simulations, as shown in Fig. [3] Below we dis¬ 
cuss how each factor in the above ratio was calculated 
from the numerical data. For recombination, the quan¬ 
tity W is calculated by averaging the quantity 
[defined below Eq. Q], in the interval [— fc#, kn]. Since 
W£ ec becomes negligible far outside of this interval, we 
define W*™ = 7 r/(fc#L) Yl k W r k ec , with 


W k 


_ reel2 

h iy 


p 0 


E|E<crM<fc<w-> ■ (ci) 


n —1 


To estimate the typical interband scattering rate from 
phonons, we evaluate W by taking an average inter¬ 
band phonon scattering rate near the resonances ±fc_R. 
The explicit form for W we used is 


W 


inter 


k R +e 


= — E 

N? ^ 


fcfl+e 

E 


w, 


k'- 


k+ 


(C2) 


fc=fc R — E k'=k R —e 


In the above, the rates W k + = W k _^( 0), defined in 
Eq. ([ 6 ]), corresponds to interband transitions from the 
upper to the lower Floquet band, through the energet¬ 
ically allowed phonon emission processes in the model 
that we studied numerically. The factor of 2 comes from 
summing over transitions at momenta around ±fcfj. Fur¬ 
thermore, in Eq. (C2), we denote by N e the number 


of fc-points corresponding to the region near k r set by 
e = The excitation density, normalized to the ther¬ 
mal density, is fitted using nonlinear least squares, to the 
form log 10 (^) =plog 10 (qf Y^S?) + ( 6 -log 10 (n t fc)) to 
obtain p = 0.49 and b = 0.95 with a standard error of 
0.001 in the region away from saturation. 


1. Scaling of phonon matrix elements 

The treatment of the phonon matrix elements follows 
along the same lines as above. Let us treat the longi¬ 
tudinal size (along the direction of the one dimensional 
system) of the phonon bath to be equal to the system size 
L. This ensures the conservation of crystal momentum, 


Appendix D: Scattering rates in the steady state 

Throughout this paper we have used the single par¬ 
ticle Floquet states to describe the steady state of the 
system. We focused on a regime in which the steady 
state approximately yields a single particle density ma¬ 
trix which is diagonal in the basis of Floquet states, (i.e., 
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the “off diagonal” correlations of the form (fka{t)fk'p), 
and higher order correlations, are negligible for k 7 ^ k! 
and/or a ^ /3). In this regime, the steady state of the 
system can be efficiently described in terms of the occu¬ 
pation of single particle Floquet states. 

Clearly, in order for the system to be in the “diagonal” 
regime described above, the lifetimes of the single particle 
Floquet states need to be much longer than the driving 
period. However, for Floquet-Bloch states near the reso¬ 
nance momenta kR, we expect a more stringent criterion 
to be necessary. To see why, recall that near the res¬ 
onance momenta, the Floquet states are approximately 
equal amplitude coherent superpositions of the original 
(non-driven) conduction and valence bands. When ini¬ 
tialized to one of these original eigenstates, the system 
oscillates between the conduction and valence band with 
a Rabi frequency equal to the Floquet gap A krt , which 
typically can be expected to be smaller than the driving 
frequency. For these oscillations (and hence the Floquet 
states) to be well resolved, and hence for our approach to 
be valid, the scattering rates between Floquet states in 
the steady state must be smaller than the Floquet gap. 

Our model includes three scattering mechanisms: 
electron-phonon interaction, radiative recombination, 
and the coupling to the Fermi reservoir. Typically in 
semiconductors, radiative recombination rates are on the 
order of Ins -1 . The Rabi frequency depends on driving 
power, and may be on the order of 0.5 ps -1 or even larger 
in experimentally accessible setups (see, e.g., Refs. (17 
and HU). Therefore, the contribution from recombina¬ 
tion to the scattering rate from a state can be signifi¬ 
cantly smaller than the Floquet gap. Phonon scatter¬ 
ing rates can be appreciably larger. An order of magni¬ 
tude for the “bare” scattering rate, (not suppressed by 
Pauli blocking), is on the order of lps -1 . However, as 
we discuss below, Pauli blocking and phase space con¬ 
siderations can significantly reduce the scattering rates 
for populated Floquet states. Finally, we study how the 
steady state evolves when we vary a coupling strength to 
a Fermi reservoir. Importantly, even when the coupling 
to the reservoir is strong enough to significantly suppress 
the densities of excited electrons and holes, it may induce 
tunneling rates to and from the reservoir which are still 
significantly smaller than the “bare” rate for scattering 
from phonons. 

When the driven system’s density is close to half filling, 
our results indicate that the steady state of the system re¬ 
sembles that of a Floquet insulator with an added density 
of excited electrons and holes. Filled states in the lower 


Floquet band with \k\ kn have suppressed scattering 
rates from phonons due to Pauli blocking. Furthermore, 
for small gy or for weak driving, the radiative recom¬ 
bination rate out of these states is strongly suppressed. 
Likewise, since these states are coupled to filled states of 
the filtered reservoir, the reservoir does not introduce any 
further scattering out of these states. Therefore, overall 
we expect negligibly small scattering rate in the lower 
Floquet band for |fc| kR. 

There are two types of momentum regions in the Bril- 
louin zone with a non-negligible occupation of electrons 
that require a more careful treatment of the scattering 
rates. Regions of the first type are those near the res¬ 
onance momenta, ±k k j, where the distributions of ex¬ 
cited electrons n e {k) = {f k+ f k +) and holes nh{k) = 
1 — () are localized. In the regime of low densities 
of electrons and holes studied in this paper, the distribu¬ 
tions n e (k ) and nh(k) are far from a degenerate Fermi gas 
due to the temperature of the phonon and Fermi reser¬ 
voirs, ksT = 0.lAfc R . Therefore, scattering of electrons 
and holes is not significantly affected by Pauli blocking in 
these regions. However, since n e (k ) and n^(fc) are highly 
localized around the minima and maxima of the Floquet 
bands, intraband scattering from acoustic phonons is sup¬ 
pressed due to reduced phase space for these processes. 
The interband phonon scattering rates are also reduced 
because of phase space arguments: recall that an electron 
in the upper Floquet band can only relax to momentum 
states around kn in the lower Floquet band; these are the 
considerations which led us to the square root behavior 
in Eq. Q. 

The second region comprises momenta between the two 
resonant momenta, |fc| < k r. This is the region in the 
Brillouin zone where the Floquet bands have an inverted 
character relative to the original valence and conduction 
bands. This, therefore, is the momentum range where 
most of the radiative recombination occurs. Since re¬ 
combination transfers electrons from full states in the 
lower Floquet band to empty states in the upper Floquet 
band, the total scattering rate of electrons in this region 
is expected to be set by the recombination rate. 

In order to compare the scattering rates in the steady 
state to the Floquet gap, we numerically evaluate the 
total scattering rate in the steady state of the system. 
Denoting by r//' the inverse lifetime of a test particle 
which we initialize in a Floquet state a at momentum fc, 
we get 


r fc°a = 5Z { W k« t 1 + N{ - £ka - + W&* a ,N{£ kla , - £ ka ) + W T k ec S kk ,6 a -d a . + } (iW) + rL [1 - D[£ ka )]. 

k'ot' 

(Dl) 


In the above, the first two terms correspond to electron-photon scattering. In the first term, the rates 
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FIG. 6. Scattering rates (red), [see Eq. (Dll] in the 

steady state of the system. The top (bottom) plot corresponds 
to the upper (lower) Floquet band. Also shown (blue) are the 
distributions n e (fc) an rih{k) in each Floquet band. In panel 
(a), we show the case of half-filling, with no reservoir coupling 
(corresponding to Fig. [3]in the main text) while in (b) we take 
log 10 (T) = 3. Other model parameters are the same as for 
the green (middle) curve in Fig. [3] Note the enhanced scale 
for the rates in the bottom plots, and the enhanced scale for 
the distributions in panel (b). 


W kc ?' = (0), defined in Eq. (|6j) , correspond to 

the energetically allowed phonon emission processes in 
the model we studied numerically. Note that W ka ^ = 0 
when £ k — < 0, due to the requirement for a nonzero 

density of states for the phonons. The second term in 
Eq. (D1) corresponds to phonon absorption, described by 
the rates W£,“, = W^,( 0) w hich vanish whe n £ k — £ k i > 
0. Furthermore, in Eq. (Dl), [see Eq. (Cl)] is the 
radiative recombination rate out of the state k, — to the 
state k , + (recall that we model recombination with “ver¬ 
tical” transitions). The results are given in Fig. [fij which 
shows rj£j normalized to the Floquet gap Ak R in both 
Floquet bands, for two representative values of the cou¬ 
pling strength to an energy filtered Fermi reservoir. The 
chemical potential of the reservoir is placed in the mid¬ 
dle of the Floquet gap. In the same figures, we plot the 
numerically obtained distributions of electrons n e (k ) and 
holes rih{k). 

Let us first examine the situation when the system is 
not coupled to the Fermi reservoir, Fig. [6jr. Consider the 
scattering rates in the lower Floquet band. For values 
of |fc| which are significantly larger than k r, the scat¬ 
tering rate vanishes; as described above this is due to 
Pauli blocking, which prohibits scattering from phonons. 
Other scattering mechanisms are absent in this momen¬ 
tum region, as explained above. In the momentum region 
\k\ < fcfj, recombination is active. For the simulations 
shown, the “bare” recombination rate (defined without 
taking into account the occupations F ka ) is taken to be 
W k ec = 3x lCT 5 A fefl , in line with experimentally accessi¬ 
ble parameter regimes. Therefore, a small nonzero l’jf?, 
set by the recombination rate, can be seen in Fig. §3- 


Tk ex ( 


Finally, for momenta |fc| ~ kjt, more significant scatter¬ 
ing rates can be observed, due to the nonzero density 
of holes and the possibility for scattering from phonons. 
However, due to phase space restrictions, the scattering 
rate is suppressed relative to its maximal possible value 
(see below), and therefore it is significantly smaller than 
the Floquet gap. 

Next, we examine the rates in the upper Floquet band. 
In this band, momentum states with |fc| ku and 
|fc| < kji are mostly unoccupied, and therefore a test 
particle initialized in these momentum states is expected 
to have more significant scattering from phonons. The 
maximal scattering rate in Fig. [6] is predominantly due 
to phonon scattering. Importantly, compared with this 
maximal scattering rate, the scattering rates at momenta 
k ~ kit are significantly suppressed due to reduced phase 
space for phonon scattering, as argued above. 

Finally, we consider the scattering rates when the sys¬ 
tem is connected to an energy filtered Fermi reservoir. 
We note that increasing the coupling strength to the 
Fermi reservoir increases the scattering rates in the upper 
Floquet band for states with |fc| < kn. This is a result 
of the significant original conduction band component in 
these Floquet states, which is coupled to predominantly 
empty reservoir states (up to thermally induced correc¬ 
tions). The rates in the lower Floquet band are only 
weakly affected, as this band is coupled to predominantly 
filled reservoir states. 


Appendix E: Particle-hole asymmetry due to an 
energy filtered lead 


In Sec. |IIID[ we studied the situation where the sys¬ 
tem is coupled to an energy filtered Fermi reservoir with 
its chemical potential /r les placed in the middle of the 
Floquet gap. Here we generically find a nonzero differ¬ 
ence between the densities of excited electrons and holes, 
An = n e — nh, despite the symmetry of the Floquet band 
structure. This highlights one of the interesting features 
of driven systems, where steady state level occupations 
depend both on the state of the bath, as well as on the 
detailed form of the system-bath coupling. 

To see how a nonzero An arises, we first note that in 
order to have /j, res in the middle of the Floquet gap, we 
aligned it with the resonant energy in the conduc¬ 
tion band. This placement manifestly breaks the particle 
hole symmetry of the system. Consider how the Floquet 
states \i/jka(t)) are coupled to the reservoir, as depicted in 
Fig- [2} Because the reservoir is energy filtered, within our 
convention for defining the quasi-energy zone the avail¬ 
able reservoir states only couple to the |^j_) harmonics 
of the Floquet states, i.e., only the rates with n = 0 in 
Eq. (101 are nonzero. Consider the situation now for 


weak driving or small gy. The n = 0 harmonics are 
then predominantly formed from conduction band com¬ 
ponents, \4>k ±) ~ | c,k). Therefore, in the lower Floquet 
band, only states with momenta —fc# < /c fcj? have an 
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Reservoir coupling, log 10 T 



FIG. 7. The offset density An. Panel (a) shows An as a func¬ 
tion of the reservoir chemical potential and the coupling 
of the Fermi reservoir log 10 T. Incompressible behavior can 
be seen when p rea « 0. (b) Vertical cuts of panel (a), show¬ 
ing An as a function of p Tea for several coupling strengths to 
the Fermi reservoir, with values indicated by the dashed lines 
in panel (a). The small slope can be attributed to activated 
behavior due to the finite temperature of the reservoir. 


appreciable |</>°_) component, and therefore only those 
are coupled to the reservoir. The situation in the up¬ 
per Floquet band is reversed: only states with \k\ > k r 


have an appreciable component, and are therefore 

coupled to the reservoir. 

From the above considerations, we see that the rates 
r° + and r°_ are not equal in the regions around the 
resonance momenta kf j. As a consequence, the following 
rates are not equal: (1) the rate for excited electrons in 
the distribution n e (k ) to tunnel to the reservoir, and (2) 
the rate for electrons in the reservoir to tunnel into empty 
states of the distribution rih{k). The difference between 
these rates leads to the non-zero value of An. 

To conclude this section, we examine how An behaves 
as fj, res is shifted away from the middle of the Floquet gap. 
In Sec. Ill D| we discussed an incompressible behavior of 
the system. Specifically, Fig. [5fi showed that n = n e + 
rih , which characterizes the number of “free carriers,” is 
unchanged by small shifts of /x res around the middle of 
the Floquet gap. Here, we complement this result by 
plotting the behavior for An in Fig. [7] The small slope 
of A n(n res ) can be attributed to an activated behavior 
due to the finite temperature of the reservoir. Note that 
|An| decreases as the coupling to the fermionic reservoir 
is increased. This is expected as |An| < |h|. 
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